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Preface 


This report is intended to acquaint engineers in their first year 
of graduate school with the fundamentals of the electromagnetic pulse. 
The mathematics is for the most part orderly starting from Maxwell's 
equations and the vector potential. In some developments, however, pro- 
cedures from advanced techniques are applied without a derivation-of 
their validity. The cited references describe these techniques. 

When I began this thesis there was concern over keeping it unclas- 
sified. I have diligently examined AFR 205-42, (U) Classification of - 
1971. It states, "Basic physics of the theoretical mechanisms by which 
NEMP is believed to be generated (is unclassified) when classified 
nuclear device design is not revealed.'"' Real and specific Jeane or 
target system characteristics or properties appear nowhere in this docu- 
ment. It is emphasized that the phenomena described here are based upon 
hypothetical, unclassified, nuclear burst parameters, e.g., the chosen 
weapon alpha (a = 1) is unclassified tecuce "Purely hypothetical values 
of a from 0.5 to 2 per shake and gamma pulse durations of a few shakes 
for an unspecified weapon may be used for purposes of illustration." 
((U) Joint AEC-DOD Nuclear Weapons Classification Guide CG-W-5, December 
1968, page 45). In my own and in my adviseor's opinion this thesis is 
clearly unclassified. 

I would like to thank Dr. Charles J. Bridgman for his guidance 
and wathusiastic interest and Dr. Gordon K. Soper for his enlightening 


comments on electromagnetic theory. 


Otho V. Kinsley 
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Abstract 


This thesis constitutes 


an introductory survey of the 


electromagnetic pulse (EMP) caused by a nuclear weapon burst. 


The survey is separated into 


(1) driving mechanism, (2) the 


electric field of the symmetric nuclear burst, (3) the elec- 


tromagnetic field from bursts whose environment is asymmetric 


due to density variations or the geomagnetic field's presence, 


(4) pulse transmission, and (5S) pulse interaction. This 


survey was collected from the unclassified literature and 


is presented here in a consistent notation and Language. 


The presentation is intended 
Standing and is not intended 
end the author has attempted 


of certain phenomena. These 


to emphasize concepts and under- 
as a design handbook. To this 
some mathematical simplification 


Simplified models are identified 


when employed and the results obtained are compared to those 


found in the literature which result from test data or from 


more exact mathematical models. 
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THE ELECTROMAGNETIC PULSE 
I. Introduction 


The primary effects of a nuclear explosion: blast, heat, 
and radiation, are accompanied by other effects caused by 
these primary effects interacting with the environment. One 
such secondary effect is the seismic disturbance in the earth, 
intense at short ranges and detectable even at very long 
ranges. Another, the subject of this thesis, is the inter- 
action of the radiation with the environment to cause the 
motion of charged particles which in turn create far reaching 
electric and magnetic fields. These electric and magnetic 
fields can be intense near the weapon and can be radiated to 
distances far beyond the range of the other effects, even | 
beyond the radiation which was their source. Because this 
radiated signal can be detected far from the burst point, 
it is important in policing the nuclear test ban treaty as 


well as an effect to be considered in atomic defense. 


Definitions 

The electromagnetic pulse (EMP) from a nuclear weapon 
may be an electric field or a combination electric and 
magnetic field which is radiated as an electro-magnetic 
Wave. in the latter case the EMP radiates 2 broadband 
Signal, i.e., one containing many frequencies, of short 
duration. Television and radar signals are broadband and 


high powered but they do not compare to the EMP in amplitude 


ry ee See aE 
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of the pulse, power output, or in bandwidth of their 
signal. 

The EMP is not the Argus effect, transient radiation 
effect on electronics (TREE's), nor the blackout effect. 
Those effects are, respectively, the input of electrons 
into earth orbit where they cause increased aurora and 
skyshine; the direct effect of the bomb's radiation on 
electronic components, €.g., transistor degradation due to 
the gamma flux; and the absorbtion of radio waves in the 
ionosphere because cf the high, bomb-induced electron 
density.? 

The effects of an electromagnetic pulse are due to 
induced currents or voltages from the interaction of the 
EMP with a physical system. These effects range fron 
ineeewenad noise in a circuit to physical destruction by 


ohmic heating, usually the former. 


Comparisons 


Although the basic physics and the magnitudes of 
lightning generated electromagnetic signals are completely 
aifferent from EMP, there sre enough similarities in their 
effects to draw some paralleis. The static from a lightning 
stroke can be transmitted around the earth by the earth's 
magnetic field. Because of the dispersion by the mediun 
through which the signal travels, the received signal sounds 
like a whistle, hence the name whistlers. Dispersion 


affects EMP wave shapes Similarly. A lightning stroke 


MiB ES ac EN 


takes place in a strong electric field which can induce 
voitages in conductors not directly affected by the primary 
stroke. The noise heard on an AM radio after a lightning 
stroke is due to induced voltages. Many of the techniques 
useful in protection against lightning strokes are also 
useful in protecting against EMP effects. 

Radar at short ranges can also exhibit many effects 
Similar to EMP. Photoflash bulbs have been fired by the 
induced currents from a radar signal. In another case the 
memory of a computer was altered by currents induced by a 
radar pulse.? 

Finally we have the radio station as an EMP simulator. 
Apecryphal stories about the early days of radic reported 
that fillings in teeth had received radio signals. In such 
a case, the filling functioned as an unintentional antenna, 
one of the hazards in protecting against the EMP. A reverse 
of the previous example is the amateur radio operator's 
dread, radio frequency interference (RFI), which is the 
spurious transmission of signals beyond the interior of the 
radio cabinets. Many of the techniques to prevent leakage 


of signals also prevents the entry of similar signals. 


History 

Although an electromagnetic pulse became a seriously 
studied problem only some time after the first nuclear 
explosion, Similar events have been studied for some time. 


One of the most widely reported examples of an induced 


electric field was Benjamin Franklin's experiment with a 
kite. Many scientists and engineers have since studied the 
effects of lightning. 

Some authors” assert that the EMP was predicted by 
Enrico Fermi before the Trinity explosion on 16 July 1945. 
An EMP was expected as evidence’ by this quote from Day of 
Trinity: 

"In the control center Sam Allison was seized by 

a sudden fear that the explosion would create a 

lighting effect and pump electrocuting volts into 

the microphone he gripped. At minus one second 

he dropped the microphone,..."4 

The same source gives us some indication of the inten- 
sity of the first nuclear EMP: 

"Scientists checking on their instruments learned 

all too quickly how devastating the blast had been. 

Only a few instruments had survived....The electro- 


magnetic storm paralysed and bent scores of expen- 
Sive gauges and measurement recorders." 


EMP Source Mechanisms 
There are three main mechanisms for generating an 
electromagnetic puise; current flow, electron-magnetic 
field interactions, and magnetic field exclusion. These 
mechanisms are described qualitatively below. Because these 
mechanisms operate differently in different environments 
there are many variations possible, but these will be 
ignored, and each mechanism will be described in a simple 
forn. 


The first source mechanism, current flow, is caused by 


energetic photons propelling electrons through inelastic 
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Scattering processes. Between the scattered electron and 
its parent ion there is an electric field. If there are 
asymetries present, unbalanced currents generate both 
electric and magnetic fields which doeether radiate a signal. 
Several sources of this asymetry are: 

a. Earth, air interface (surface burst} 

b. Case asymetries (case signal) 

c. Exponential atmosphere (medium altitude burst) 

d. Space, atmosphere interface (high altitude burst) 
For all these asymetries gamma rays are more important than 
are X-rays except for a height of burst cbove * 100 kilo- 
meters. 

The second major mechanism, electron-magnetic field 
interaction, operates in the presence of the earth's mag- 
netic field or the field generated by currents due to the 
explosion, In a dense medium the interaction of electrons 
with a magnetic field is a secondary effect but does occur. 
In the thin atmosphere of near space the interaction of 
electrons with the earth's magnetic field is very important. 
There are two regions of interaction: around 20-30 kilo- 
meters where gammas and neutrons are absorbed and around 
100-110 kilometers where X-rays are absorbed. In these 
regions the earth's magnetic field causes the electrons to 
move in circular paths, and these accelerated electrons 
radiate. The radiated signal is energetic and has a low 
characteristic frequency which is important because low 


frequencies propagate for long ranges. 
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The third mechanism is the exclusion of the magnetic 
field by the expanding ionized sphere. The compression and 
relaxation of the magnetic field lines radiate a signal, 
but it is minor compared to the other sources. 

In passing we note that the regions of interest and 
mechanisms can all te considered axially symmetric if local 
air density and topographic variations are ignored. Because 
the deviations from axial symmetry cause only second order 


effects, axial symmetry will be assumed in order to simplify 


the mathematical description of the problem. We will discuss 


this point more fully later. 


Transmission of the EMP 

The radiated EMP is like any other radiated signal in 
that its propagation can be divided into a near field, an 
intermediate region, and the far field. The near field may 
contain a non-cadiated component and components which are 
attenuated rapidly. In the far field only a dipole compon- 
ent is significant.. The intermediate region is a combina- 
tion of the two. 

The electromagnetic puise can be thought of as the 
summation of a series of trigonometric functions of differ- 
ent frequencies, i.e., a Fourier series. Because a signai's 
velocity of propagation is usually a function of frequency 
the shape of the EMP will change with distance. 

The signal may be attenuated by the medium. For 


example, in earth, representative figures show that a 
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Signal will be decreased by lfe in about 100 feet. In air 
the attenuation is selective, being more rapid for higher 
frequencies. In the ionosphere the attenuation is dependent 
upon the electron density which may have been affected by 
the nuclear burst. Because an analytical solution for the 
electric or magnetic fields is virtually impossible, many 
studies involve numerical solutions to the field differen- 


tial equations for different situations. 


Reception of the EMP 

The final area of interest is the reception of the EMP. 
by systems which might be damaged by it. The EMP can enter 
a susceptible system by induced currents or by conduction 
currents from magnetic and electric fields respectively. 
The effects are proportional to such things as area, resis- 
tivity, and a coupling coefficient. Some examples in 


Chapter VII will illustrate reception, 


Overview 

The next chapter reviews Maxwell's equations because 
they play a vital role in understanding electromagnetic field 
relations. The reader who is well grounded in electromagnetic 
theory may skip this chapter. The succeeding chapters cover 
the source of the currents, the symmetric burst, asymmetric 
bursts, transmission of signals and reception of signals. 
The asymmetric burst is discussed under three major subdivi- 
sions: asymmetric current flow, magnetic field-electron 


interaction, and magnetic field exclusion. 
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Each chapter is introducea and summarized. The 
equations are numbered serially. The appendix details the 
mathematics behind several processes. A fold-out, page 191, 
is added to allow Maxwell's equations, conversion factors, 


and some basic data to be available at all times. 
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Il. Maxwell's Equations 


Maxwell's equations describe the behavior of electrs- 
magnetic phenomenon in almost all their details (any excep- 
tions appear as a result of quantum effects and do not 
affect the understanding of EMP}. ‘he four Maxwell equations 
Will be introduced through mathematical manipulation and 
physical arguments. After being introduced they will be 
used as if they were the basis of electromagnetic phenomenon 
without regard to their development. 

Actually Maxwell originated only one of the four 
equations, but with it he unified the work of the persons 
whose names are associated with the other three: Gauss, 


Ampere, and Faraday. We will consider these in order. 


Gauss' Law 

Before considering Gauss‘ law, it is necessary to state 
the work of C. A. Coulomb. Coulomb formalized the experi- 
mental relation between two electric charges as a force 
proportionai to the magnitude of each charge and inversely 
proportional to their separation distance squared. 


Coulomb's law is 


Pree 
F YT, VF 


(newtons) (1) 


where 9, 5 es are the charges in covulombs; Ff is the line 


between 9: and gz of length * in meters, and €, is the 


Law 


Gauss' 


10 


permittivity of free space (here a proportionality constant) 
in units of farads per meter. In otder to consider the 
effects due to an isolated charge te introduce the electric 


field vector, fal » the force on a unit positive charge: 


Ea lm = BEWESES gu VOUtS,. 19) 


Now Gauss' law relates the total charge within a volume 
to the total electric field vector, E » penetrating the 
bounding surface. Consider a single charge 4 which is 
enclosed inside of an arbitrary volume V as shown in Fig. 1. 
The total surface area of this volume is designated S - By 
Coulomb's law, the electric field at any point on the surface 
is 


ae 9 r 
Fe ore (3) 


where [ depends upon the location of the charge ys and upon 
the point of the surface which is under consideration, 
Uniess the voiume is a sphere with 4 at the center, F 


changes from one surface point to another. Now let us 


consider only the field components normal to the surface 


— vi r:-n 


476, 73 


> 


(4) 


and sum this equality over the entire surface of the 


arbitrary volume: 
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PE Ads =O 2 (5) 
s A 


yIre, 1? 


— A 
The factor 3 /4T € commutes andY-NM may be replaced by 
3 & A . 
rcos O where G@ is the angle between [ and n , - 


AT & re 
But, as can be seen from Fig. 1, cos 9 ds = da where dA 


is perpendicular to I ana 9A/,* is a differential increment 


of solid angle,df2 . The summation over all solid angle is 


41 , so 
7 ° 
‘S ° 


By superposition, we can extend this argument to a number of 
charges, a > or to a volume distributed charge, CF) 5 


coulombs/m°, in which case 


E-ntds -[ 72 dives (8) 


s 


which is Gauss' law. 
Gauss’ law can also be expressed in differential forn. 


A theorem due to Gauss himself states that 
gma ds -[ vA en (9) 
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Thus the ieft hand side of Eq (8) becomes / V~* dvew 


and we have 


(v-2- 22) dum = (10) 


But since the volume is arbitrary, Eq (10) can be zero only 
if the integrand is zero, which yields Gauss' law in differ- 


ential form 


wee ee! (12) 


Eo) 


which is also the first of Maxwell's equations.& 


The electric field vector, Ee ,» and © may be combined 


to yield 


VD =AlF) - (12) 
where D is known as the displacement vector. We shall 


return to the divergence of D during the discussion of the 


physical meaning of Maxwell's equations. 


Ampere's Circuital Law 

Another of Maxwell's equations first appeared as 
Ampere's circuital law which states that the magnetic field 
around a closed loop eauals the total current through a 
surface bounded by the loop. We shall derive an integral 
Statement of this law using the force between current 


— 


elements and a magnetic field vector B which has much the 
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Fig. 2. Ampere's Law. 
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same fundamental role in a magnetic field as EF does for the 
electric field. 

Ampere, around 1820, studied the forces between rigid 
conductors for various values of the current they carried.? 


He arrived at an inverse square relation of the form 


F = ==. z ap a ae R) (13) 
Mes lal a R? 


constant of proportionality (magnetic suscepti- 


= 
i 
o 
+} 
ti 
> 
o 
HT 


bility 47x 107 henry/meter). 


N 
Ky 


current in circuit (amperes). 


Qn 
“1 
io 
bal | 

f 


= incremental length of circuit with conventional 
orientation, i.e., bounded surface is on left 
of direction of travel. 
Renee lee oa te eek 
Reyer = vector from dr to dr of magnitude (meters). 
f « - 
and C and C are any two closed circuits such as shown in 
s 


Fig. 2. The increment of force that all of circuit C 


exerts on an increment of circuit C is 


inj. LD dex Ge«R) 
YT fe" rad R° 


Just as in the electrostatic case where E=F/Q , let 
us define a magnetic-force-per-element-of-current as a field 
quantity, B , which is called magnetic induction (in 


webers /meter*) > 


ee BF) = d Fron (15) 
an SS 
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or by Eq (14) 


= ee ae 7 
Bu = Ant (16) 
41 fs E 
fae. OR 


Compare with Ew) = rere _ 


where the primed coordinates refer to the source of the 
field in contrast to the unprimed coordinates which are 
field points. | 

By using vector identities and the properties of the 
del operator, V » we shall derive three useful relations: 
continuity of the B vector; a vector potential A 3 and 
Ampere's circuital law. Since WR) = Rife the expression 
for B becomes , 


= r fos cL : 
B= l i dr’x VER (17) 


By the vector identity dex v{¢] = Vx di’ eats vx dF?’ 


Eq (17) becomes 


¢ 


but the last term is zero because the del operator is over 


fou 


the unprimed (or field) coordinates 
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ee 


yy | T's x dr’ 
or Be Tie pee (19) 
Cc 


Further, the del operator comm»tes since it still does not 


affect the field-source terms, 
thus B=Vx = = (20) 
47 f, “RR | 


By taking the divergence of this expression we see that, 


- } eo | 
2B Vv ate 4 d7 
V:B=V-V: | FIA f = | (21) 


but because of the vector identity V-Vx{ )=0 


V°-B=0 (22) 


Equation (22) is a mathematical statement of the continuity 
of the B vector and is another of Maxweli's laws. 

The bracketed function in Eq (20) is designated the 
vector potential, A - Now B=V x A in all cases. From 


vector analysis 
Vx( ) =-V"( + V(0-0)) (23) 


which, for the magnetic induction has a zero final term; 


therefore 


V*B=-V’A (24) 
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Carrying out the indicated Laplacian operation on A gives 


VxB = = (25) 

or 
VxH= F (26) 
5 2 


where J is the source-current density in amperes-meter™ 


and 


H- pe B (27) 


where H the magnetic field in amperes/meter. 


Summing the current density through a surface S bounded 


by a curve C leads to Ampere's circuital law. Integrating 


Eq (26) 


Fd ed vxH-d3 (28) 
AS 5 . 


and by Stoke'ts theorem 


1 


vxH-d5 = H-d2 (29) 
> 
therefore 
_~ j ri 
L - H -a f (30) 
where I = [F-d3 is the total current through the surface. 
Is : 


Equation (30) is Ampere'’s circuital law, the steady state 
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form cf one of the Maxwell equations. The general form 


will be shown after noting two time dependent relations. 


Faraday's Law 

Faraday's law, the third of the Maxwell equations, was 
developed about 1852 by M. Faraday and, independently by 
J. Henry, but by right of first publication Faraday's name 


10 


is commonly associated with the effect. These men experi- 


mentaily found that the voltage induced in a circuit was 


proportionai to the rate of change of magnetic flux 
=- (31) 


where Ui is the total magnetic flux in webers through the 


surface of the circuit: 
w= | B-ds (32) 


and 


‘dl (33) 


< 
" 
mi 


=) 


where V is in volts. From Eqs (31), (32), and (33) we have 


the integral form of Faraday's law)! 


Eide Ss g_[B-a5 (34) 
dt 
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For the differential or point formulation we tra ‘sform 
the line integral to a surface integral by Stoke's theorem 


and have 


cr 


jvxé-ds 2 ua -d5 (35) 
+ 


where Se Ss) or 
[Vx +8 |-d5 =0 (36) 


Because the surface S is arbitrary, the integrand must be 


zero and 
VxE + 28 210 (37) 


which is the differential form of Faraday's law and another 


of Maxwell's equations. 


Conservation of Chares 

Before we examine a gedanken experiment which illus- 
trates Maxwell's contribution, Jet us consider conservation 
of charge or the continuity equation. In a volume V, the 
net change of charge is Ye which must equal the flow into 


or out of the surface, given by V-J or 


a2 = -V- (38) 


dt 


U 
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The minus sign is present to account for the outward normal 


to the surface. 


Maxwell's Contribution 

A discrepancy which J. C. Maxwell's contribution 
removes can best be seen by applying Ampere's circuital law 
to surfaces I and II in Fig. 3. The current which flows 
into the capacitor for a time has a magnetic field associated 
with it. The magnetic field found by considering surface I 
is H= lvg 3 but considering surface II H=0O since there 
is no current through that surface. In 1861 Maxwell], proposed 
a current, GD (called displacement current}, so that 


Eq (26) becomes 


a 
| 


VxH ae 


Qa 
ck 


which experiment soon proved to be correct. This displace- 


ment current remeves the illustrated discrepancy.+3 


Maxwell's Equations 


A summary of our results shows the Maxwell equations: 


I VAD) — 8) =O (39) 


Il V°B —O (49) 
TEI VXxE + 32 pclae oS (41) 
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Fig. 3. Maxwell's Contribution. The current loon 
penetrates two surfaces. Surface I is penetrated by 
the wire of the circuit. Surface II is not penetrated 
by the wire but instead passes through the capacitor's 
plates. Surface I and II have a common boundary. 
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ete Mi 


IV — Vx - 3 ar =O (42) 


and two additional equations: H = fro B and 
D-2G.E- 2 

Note the similarity between the equations and that E 
and BD hold places similar to B and H respectively; also 
note that €, and pe" have similar functions. This can aid 
in transferring mathematical ideas from one type of field 
to the other. 

In Eq II we see that there is no magnetic pole distri- 
bution comparable to the charge distribution in Eq I. That 
is a statement that no isolated N or S magnetic poles exist 
or every point in space is magnetically "neutral". This has 
its parallel in Eq III where we see no magnetic current, 


JN-poles or d5- poles is evident. 


Word Picture 

A word picture of the physical meaning of Maxwell's 
equations is possible. 

Equation 1,V-D -A( a © » implies that the D 
vectors start on positive charges and end cn negative charges; 
in the presence of material such that D - c. E » the 
ends of the F vectors can be thought of as associated with 
induced charges, and D vectors continue to be attached to 
real charges.?4 (See Fig. 4). 

Equation 11,V° B= QO , states that B vectors form 


continuous loops. The distinction between B and H is 
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Fig. 4. Tllustration of Continuity of Field Lines. 
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similar to that between E and D but with H being created 


by induced or amperian currents and B due to true currents. 
H vectors start and stop on isolated poles which can be 
imagined for their mathematical convenience. 

Equation 111, PE-dl = -}, [Beds (integral form), 
shows that a changing magnetic field creates an electric 
field. 

Equation wv, g Hdl = 2 [6 d5 + [7-45 {integral 

dts As 
form), shows that a changing electric field creates a 


magnetic field and. so does a steady current. 


Boundary Conditions 
in-general the boundary conditions on the field vectors 


are that the normal component of B is continuous, the 

tangential component of E is continuous, the normal con- 
ponent of D changes by the value of any boundary surface 
charge density (@5 ) and the tangential component of H 


changes by the value of the surface current density az y 35 


In equations we have 


B.* By =O (43) 
Ey - Ey = 6) (44) 
Da: Dep ee (45) 


Hy -Hy= Jy (perfect conductor) (46) 
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Two additional relations are heipful in dealing with 


electromagnetic preblems: the Lorentz force and Ohm's law. 


Lorentz Force Equation 

The Lorentz force equation is easily derived from the 
theory of special relativity. (The equation's second tern 
is a relativistic correction to the force qe 5.18 Here 


we shall only present the equation and suggest that it is 


true; it is 


F=qE+qv«B (47) 


- 
~ 


where V is the velocity of the charge Q relative to the 
field B. In the absence of an electric field Ee » the 


equation is (for a differential of charge dq ): 


GF 


= dq Vv xB (48a) 
= dq drFxB (48b) 
di 
=dq drxB (48) 
dt 
= L dFxB (48d) 


which is the vector form of the defining equation for B. 


Thus we have indicated the validity of the Lorentz equation. 
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Vesthm.s aborted 


Ohn's Law 
Finaily we have Ohm's law, V =IR or as we shall 


use it 


J=ocfk (49) 


Conclusion 

Maxwell's equations and the constitutive equations 
presented here are slightly restricted in that some linear 
relations (specifically D= 6& and H = u,'B) were used 
whereas reality is usually non-linear. These restrictions 
should not affect the material as presented in this thesis. 

In general, the vector and scalar potentials are easier 
to deal with than are the E and B vector relations; the 
point forms of Maxwell's equations are more useful than the 
integral forms; and adroit use of physical insight, well 


chosen boundaries, and appropriate simplifications are 


needed to solve Maxwell's equations. 
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EII. Current Source 


Overview 


The electric and magnetic fields of the EMP are gener- 
ated by currents near a nuclear burst. These currents are 
the motion of electrons or ions outside a sphere which con- 
tains a nearly full ionized plasma in which there is no. net 
current flow (This isothermal burst sphere is the fireball 
of early times up until hydrodynamic separation*). The 
initial motion is that of electrons moved by energetic pho-. 
tons by Compton scattering hence the name Compton current. 

After a photon-electron interaction, momentum must be 
conserved and since the photons move radially outward, so 
too will most of the electrons. This movement of electrons 
is a current, the asymmetries of which dictate the electro- 
Magnetic field to expect. For example, an external magnetic 
field (such as the earth's geomagnetic field) causes the 
electron to turn generating currents which subsequently 
radiate an EMP. Bigure’s summarizes the process. 

The behavior of the currents generated depends largely 
upon the behavior of the photon sources. In this chapter we 
Will consider photon sources (gamma sources, X-ray sources and 
nmeutrons as a source}; the photon attenuation; and finally the 


relation of photon flux to current. 


*Glasstone, Samuel. The Effects of Nuclear Weapons, AFP 
136-1-3. Washington: Department of Defense, April 1962. 
(Page 71) 
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Fig. 6. 


Time Dependence of Bomb X-Ray or Gamma Output.18 
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Photon Production 

Figure 5 shows that several photon sources are present. 
We will discuss the time and magnitude behavior of gamma 
rays, X-rays, and neutrons in the following sections. 
Simplified models of the sources will also be discussed. 

Gamma Rays. The gamma source can be modeled by a point 
source whose intensity grows exponentially to a peak, then 
decays exponentially at a slower rate for a time, then is 
constant, and finally drops to zero. This is illustrated 
in Fig.6. For some uses the central source is modeled by a 
sphere of finite radius. In this thesis we will assume the 
gamma source to be 

GN. 2. eI) (50) 
t E 

where Y is the yield of the explosion (in MeV), F is the 
average photon energy (MeV), € is the fraction of the yield 
released as gamnas, f(t) is the time variation of the source, 
and N is the number of gamma rays. 19 This time variation 
is normalized so that 


oo 


wate: 4 (51) 


— oo 


The assumed time variation is (on the basis of Fig. 6) 


f(t) =O , te (-00, O) (52) 
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f(tj= 0 té (T,,0c°) (53) 


f(t) = Aext te LOT] (54) 
fi) = Bert-D sve(T, G3 (55) 
f(t)=C rté (GI (56) 


where we chose 1, ~ &:1,~30T= 7% ; 1, ~300T, 


= 900% and 104 = o . Through the boundary 
conditions that f(7,) = (7) we have 
B= Aexh (87) 
Cspe eV (58a) 
=Ae a MAN) (58b) 


By normalization we have AS “Yo95 

For Example 1 let us find the gamma source function of 
a one megaton explosion which emits 0.3 percent of its yield 
as 1.5 MeV gamma rays. We have 


(59a) 


aN 6 3 (ensieneee NENG 003)10% € ot 
1.5 (MeV) . 225 
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or 


die 2 Bix 1071 ent gammas AclOT] (59b) 


——$_—_— 


t second 


where we have chosen a hypothetical alpha of 10*® second7?. 


The Spectrum used in the preceding example was a mono- 


energetic photon source. Using the Maienschein prompt gamma 


spectrum, R. E. LeLevier found that at sea level the main 
contribution to the Compton current was due to 1.5 MeV 
gammas out to 500 meters and due to 2.5 MeV gammas at 1000 
meters (Fig. 7).70 Therefore we shall use 1.5 MeV as the 
mean photon energy throughout this thesis. 

This mean energy is due to a combination of prompt 


fission gammas and inelastic scatter gammas. The spectrum 


of fission gamma rays has a mean energy of ~~ 1 MeV; al MeV 


gamma implies an average forward-scatter energy of 0.4 MeV 


for the Compton electrons. Because fission and fusion neu- 


trons have energies which exceed the threshold energy for 
inelastic scatter with case materials, a second, distinct 


spectrum of gammas is produced which has a mean energy of 


4-5 MeV. The average electron forward energy from 4 MeV 


gammas would be ~/2.S5 MeV. The ion pair production due to 


a 1.5 MeV photon is “2 x 10% ion pairs per photon (see 


Eq (99a)) which compares well with 3 x 10° cited by Karzas 


and Latter.?} 


‘ 


The gamma spectrum grows harder farther from the source 


because the higher absorption at low energies removes the 


low cnergy photons. The spectrum grows harder with time 
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Alnterpolated 
points 


Compton current/unit energy. interval (arbitrary units) 


Energy (Mev) 


Fig. 7. Compton current per unit energy interval 
for the Maienschein spectrum. ? 


34 


Sears Sone 


GNE/PH/71-4 


Since the average energy of the nitrogen capture gammas is 
~5 Mev, 2? much higher than the assumed initial average of 
1.5 MeV, 

X-Rays. The extremely energetic particles and contained 
heat energy within the interior of the fireball quickly 
equilibriate and are then the equivalent of a blackbody. 


Karzas and Latter> 


state that the X-ray output can te 
modeled by a 1000 eV (1 KeV) blackbody; the mean energy of 
a 1 KeV blackbody is 3.8 KeV. About 80 percent of the yield 
of fission weapons is emitted as X-rays. We shal} assume a 
1 KeV X-ray spectrum and assume that the X-rays follow the 
same time and space variations as gamma rays, Eq (50). 
Neutrons. Prompt fission gammas are the major source 
of early gammas, but at later times neutron inelastic scat- 
tering in the air, and neutron inelastic scattering and cap- 
ture in the ground produce most of the pangs. Still later, 
the isomeric decay of fission fragments, neutron capture in 
the air, and, finally, gammas following beta decay of fission 


products dominate the photon source. 24 


In space the proc- 
esses depending on the presence of air must be neglected. 
Neutrons interact in several other ways, e.g., (n,d)3 (n,p)}3 
and elastiec scattering which can produce an energetic ion. 
These may require study in a detailed calculation of the 


EMP.2°> 


All of the interactions listed would produce intense 
ionization and affect the conductivity of the region in which 


they occur. 
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GAMMA OBSERVER “*¥ 
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a 


Fig. 8. Stationary Neutron Source Model. The observer 
at P can approximate the gamma source as a point until 
the neutron wave front is comparable in size with [. 
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A typical inelastic scatter reaction is the nitrogen 


reaction 


which may occur when the incident neutron energy exceeds 
the threshold energy. The cross section for this reaction 
is about 0.5 barns if the neutron energy is S 7 Mev and 
4-5 MeV photons are emitted. 7° 
The gammas produced by neutrons account for the late 
time behavior of the source. This neutron source can be 
modeled or approximated by the "rattle" or stationary source 
model which assumes that all neutron induced gammas originate 
from the burst point itself. {See Fig. 8.} In this model, , 
neutrons rattle about near the burst point being absorbed 
according to their velocity V and inelastic scatter cross 


section Wilareue? This results in a gamma ray production 


‘ rate of 


Q. 


x = N(t), v (n,n*) (60) 


where N(t), , is the number of neutrons present at time 35 
oO 


and is 


Nt), = No e Lvt (61a) 


iY (61b) 
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width = v At 


Fig. 9. Spherical Shell Model for Monoenergetic 
Neutron Pulse. Current and ionization rates at P can 
be calculated by the use of spherical harmonics. The 
neutron shell is modeled by a square pulse. 
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or the gamma production rate is proportional to ar 27 
This approximation is good while the radius of the 

expanding shell of neutrons is small compared to the dis- 

tance to an external observer. A 14 MeV neutron has a 


7 meters-sec™} 


velocity of 1.52 x 10 or ~~ 1/6 the velocity 
of light so we see that the gamma flux does precede the 
neutrons for early times. 

A model valid at later times is the shell model which 
assumes that the neutrons are in a shell whose radius 
increases at the velocity of tne neutrons. Shaeffer and 
Hale describe the shell model which is illustrated in Fig. 9. 
The current and ionization rates will peak when the neutron 
shell arrives at the observer point due to the decreased 
effect of the Vp2 attentuation of the gammas emitted from 


the shel1.29 


Shaeffer uses this model for a spherical shell 
of monoenergetic neutrons; Hale? compares this model with 


more accurate calculations. 


Photon Attenuation 
The virgin or unreacted flux from a point source is 


given by 


D (es IND) G(F) (62) 


or 


(dr 
Ort Y= anyed = 


2 ; r 
di. 4%r ieee 
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Meters“/Kg Air 


Photoelectric 
Absorption 


Total 
Absorption 


Compton 
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Pair y 
Production—~*/ 
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/ 


Fig. 10. Mass Attenuation Coefficients for 
Photons in Air. The curve marked “total 
absorption” is (74) =) +) +(*p) 

where od ,7 andk are the corresponding linear 
coefficients for Compton absorption, photoelectric 
absorption, and pair production; “ is the air 
density. For example at 1.5 MeV the Compton 
absorption is 2.4 x 107 meter*/kg of air; at 

0° C and 760 mm Hg pressure, the density of air 

is f = 1.293 kg/m3.51 
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where [. is radial position eekesy and [L{T) is a 

cross section for photon removal (meter~!),32 The assump- 
tions that allow this function to be used for all purposes 
include time independent cross sections, no scattered photons; 
and no time retardation, 

These assumptions are valid for a rapidly increasing 
source rate; one much larger than the rate of change of 
electron density due to scattered photons. If the source 
rate is rapidly decreasing the Compton current and energy 
deposition rate are mainly determined by the scattered 
photons, i.e., the photons associated with the build-up 
factor. If the time variation of the source is slow rela- 
tive to the time-spreading induced by the multiple scattering 
then both virgin and scattered photon reactions must be 
considered.?°> 

The cross section for photon removel, fee » is the sum 
of the cross sections for three processes: photoelectric 
effect, Compton effect, and pair production. Each is domi- 
nant in a different energy region as shown by the graph of 
their cross section, Fig. 10. X-rays are mainly attenuated 
by the photoelectric effect and gamma rays are most heavily 
attenuated by the Compton process (for the energies we are 
concerned with}. The forward scatter of photoelectrons is 
less pronounced than the forward scatter of Compton electrons 
because the Compton effect does not involve the nucleus and 
the high gamma momentum is partially transferred to the 


electron. 


4] 


At high energies ( > 1.02 MeV) pair productior 
dominates. The two particles produced have a radial compo- 
nent to their travel; however they do not contribute to an 
electric field since they are neutral overall, i.e., an 
electron and a positron cancel in charge. We shall not 
consider pair production further. 

Photoelectric Effect. The photoelectric effect is due 


to the force upon an electron due to the electric field, E 


> 


of the incident photon. Because Eis normal to the incident 


photon's direction of travel the predominant scatter at low 
energies 1S normal to the photon momentum vector. Because 
all the energy (less the binding energy) of the photon is 
absorbed by the electron,the electron's momentum is greater 
than the photon’s. To conserve momentum the atom recoils.>4 
Because a bound electron is in an energy well, the 
incident photon's energy must exceed the well's energy 
deficiency in order to free the electron. This implies that 
the cross section for photoelectric effect will take sudden 
jumps as the photon energy increases enough to free more 
tightly vound electrons from an atom's inner shells. The 
increase in cross section is due to the increase in the 
number of available interaction particles. Of all the major 
air constituents, oxygen has the largest absorption edge at 
0.532 KeV, well below the average energy (3.8 KeV) emitted 


by a 1.0 Kev blackbody.°°> 


Pere 


Overall, however, the photoelectric cross section 
decreases as photon energy increases. The decrease can be 
thought of as due to the decrease in the time available for 


the E vector to accelerate the electron (the frequency 
increases as the photon's energy increases so the E vector 
direction changes more rapidly). 

The angular distribution of the photoelectron becomes 
more forward biased as the photon energy increases, but the 
biasing is not really strong until the photon energy is in 


the 10's of KeV. An expression relating the electron for- 


ward momentum and incident photon momentum is 


\% Average Photoelectron 
Re 442 cos @>= Forward Momentum (64) 
ame Photon Momentum 


where R is ~w 1.5 - 1.6 for photons with energies from 

0.5 MeY to zero2é Karzas and Latter?’ assign the value 

AG oh to the mean cosine of the scatter angle, 7, 
fer the photoelectric effect. For photoelectrons with 1] to 


a few KeV energy they state that 7}? 0.05 to 0.1. 


Compton Effect 

In our discussion of the Compton effect we shall ignore 
scattered photons and concern ourselves oniy with the 
electron dislodged by a virgin photon traveling radially 
outward from a central source. The scattered photons are 
ignored because their energy is about half the incident 


photon's energy and because their direction is somewhat 
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Fig. 11. Compton Scattering. 


- Electron Energy (MeV) 
fh 


/ 


e) 30 50 90 
Angle (degrees) 
Fig. 12. Compton Electron Energy as a 


Function of Scatter Angle, Eq (68), for 
Several Incident Gamma Energies. 
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random. Because the electrons scattered by these secondary 
photons would be directed even more randomly, their main 
contribution to EMP would be to increase the ionization. 
Since E = S/o and the conductivity G@” is a function of 
the ionization we see that E would decrease if O increases. 
Using a buildup factor would result in too large an estimate 
because a buildup factor has no direction associated with it 
whereas the direction is needed to calculate the current 
preduced. We propose to look at the maximum threat there- 
fore secondary photons will be ignored. 

In Fig. li, h 2% is an incident gamma ray and hv is 
the scattered photon. The electron wrnich was initially at 
rest at © is given a momentum P - Because this is a two 
body problem it can be worked in a plane. Conservation 
laws and the relativistic expression for kinetic energy 


allow us to find an expression for the electron kinetic 


energy as a function of scattering angle and photon energy: 


Az = 

= = Ne + P (Conservation of Momentum) (65) 
hv, = Ne « JT (Conservation of Energy) (66) 
Pc = Eh 2mec (Relativistic Energy) (67) 
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where 
P Electron momentum 
C Velocity of light 
MN Planck's constant 
hr, Incident Gamma Energy 
NV Scattered Photon Energy 
T Electron Kinetic Energy 
MM, Electron rest mass 
Solving Eq (65), (66) and (67) for [| and, 


after Evans, letting & = Reyes , we have 


T= CNBEX cos“® (68) 


 (qaocye - <cos*O 


which is plotted in Fig. 12 for several incident photon 
energies. 

The cross section for an electron to be scattered into 
a cone with walls at W and o+d9 | is plotted in Fig. 13. 
Note the strong forward scatter. If this odes section is 
altered by an energy fraction, Yow, » the plot is even 
more forward biased, Fig. 14,38 

The average energy of the forward scattered electrons 


was computed using Be 
a 


TD) RCP) dp 
<T> = az 
pd) db 


(69) 
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Fig. 13. Number-vs.-angle Distribution of conpeon 
Electrons for Several Incident Gamma Energies. a 
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Fig. 14. Compton Effect Energy Distribution. Differential 


Cross Section per unit Angle for Electron Energy Scattered 
in the Direction b 41 
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Fig. 15. Average Electron Energy and Forward Energy 
as a Function of Incident Photon Energy. 
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where 
CT? The average electron energy 
TD) The scattered electron energy as a function of 
electron scatter angle 
(Pp) The cross section for scatter into an angle 
p> Prag 
The average energy flow parallel to the direction of the 


incident gamma was found by 


y 
T(P)cosH p(d) dh 


[ucores 


where Tp)cosh is the component of the electron 


Theo? = (70) 


energy in the forward direction, I. 


Ewp - The results of 


these two calculations are shown in Fig. 15. The forward 
energy of the electron created by the “standard" 1.5 MeV 


gamma is 0.7 MeV. 


The Compton Current 

Cross Sections. Because the relations developed in 
this section for use throughout this thesis invoive the use 
of cross sections and mean-free-paths, a few paragraphs will 
be devoted to discussing them. Many EMP models concern air 
and its variation with height, e.g., the earth's surface may 
be represented by air a thousand times as dense as air at 


stp. 42 
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An effective microscopic cross section a LL 
(barns/particie), can be assigned to air because it has a 
fairly consistent composition below 100 kilometers. The 


macroscopic cross section,* {L (meter-1), is 


fh = NCA) oe (71) 


where N(N) is the number of particles per unit volume 
(= P4 n Na G awonyi particles/m°> and NM is height. The. 


ratio of two cross sections does not vary with altitude: 
phe 2 NCh)abé 2 oe (72) 
- type i= 
bj Nh) By abts 


Because of the reciprocal of total cross section is photon 
mean-free-path Ay the ratio of mean-free-paths is constant 


also: 


Ba = Wy oe 
By ”, Ki 


The macroscopic cross section at a particular height 


(73) 


iT 


can be related to the cross section at a base height through 


the relative density, P/p, 


*The cross section, by , in meters-~! is related to the 
cross sections of Fig. 10 through the density, i.e., 


plot!) = p(m?-kgt) p(kg-m3) 
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PNA 
Bi GAW ~ aft (74a) 
a 
= ~~ Hy, (74b) 


The mean-free-paths are related by 


In the discussion of the Compton process the total 
cross section and the Compton cross section were equated to 
simplify the discussion. In fact the photon flux attenua- 
tion is proportional to the total cross section, fed » but 
the Compton electron production is proportional to the 
Compten cross section, © . The difference is small for 
1.5 MeV photons. 

Another anomaly occurs in the treatment of electron 
cross section and range. Because an electron can Scatter 
through a large angle while stopping, its mean free path is 
usually much greater than its mean range (termed mean-free- 
range (MFR,)). For example, the electron distribution within 
a material one mean free path thick which is bombarded from 
one side by an electron beam is greatest about half way 
through the material. For this illustration the MFR, is 


about half the mean free path, Reo? 
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Fig. 16. Electron Fiux Production. Electrons produced 
within one MFR, of the right face by a constant gamma 
flux from the left will penetrate the right face. Since 
the right face has unit area this electron movement is 
just the electron flux. 
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The Compton Current. The net forward scatter of elec- 
trons is a flow of current which is called the Compton 
current, J - The rate of generation of free electrons obeys 
the general reaction rate formula, fiux times macroscopic 


cross section: 
Electron Production Rate = Pir, t) oO (76) 


where Py is the instantaneous photon flux at the point 
and O is the Compton cross section. 

A pictoral description of the production of a Compton 
current supposes a plane wave of photons incident from the 
left onto a series of rectangular boxes of unit cross 
section and stacked deep enough so that any electrons formed 
in the volume will escape out the right end, i.e., the stack 
is one electron mean-free-range deep (Fig. 16). Assuming 
that the photon flux does not change over this distance we 


have 


Paz PyMFR, (77) 


i.e., electron flux is photon-flux times path-length times 
_ cross-section-per-unit-path-length. Now current density is 


electron flux times electron charge 


a 


Pa 4, (78a) 


or 


dyMFR./ €A, (78b) 
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J. = h—£e4 (78c) 
Cc , >» r 

t 
where we have replaced the Compton cross section O by the 
total cross section LL (= Vy which assumes that Compton 
Scattering is the only reaction taking place. 

If we define an electron cross section x = 1/MFR, 

we can re-write Eq (77) as Pade = Ly Py - Transposing 
and multiplying by © produces an expression which is the 


rate of change of the charge density: 


oO 


SH = e(beLe ~My) = 0 (79) 


This is proportional to the divergence of Je $o our develop- 
ment is equivalent to specifying V’ Je=0 from Eq (38). 

_ Example 2. This example illustrates the magnitude of 
the Compton current. The sample nuclear burst has, fron 


Eq (59), a maximum gamma flux rate of 


«xt 
= = 23x107e tet (80a) 
- 46x10°° 1.5 woe ganmas (80b) 


“ 


Let us find the Compton current at a point one kilometer 
horizontally from a burst at a height of ten kilometers. The 


data concerning 10 km altitude is found in Table 3, page 96 


$4 


With srurgr ame 


and is [Ly = i.11 km™1, Ay = 900 m and MFR, = 9 m. 


flux at the point is 


E 


$= SSar) 


dt 


= 46x1082 2 
ATR 


z 1.2 x10e° 1.5 MeV gammas 


mé- sec 


From Eq (78) the electron flux is 


MER 
Pe 7 oF te 


25 9 
= 2xl0> == 
SOO 
- 12 x1023 electrons 
m*-sec 
The Compton current is 
Jo = Pee 


= (1.2 x10°" 16x10") 


55 


(62) 


(81a) 


(8lb)- 


(78c) 


(82a) 


(82b) 


(78a) 


(83c) 
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J =2xi10% amps (83) 
Across a surface about the size of the lead in a wooden. 


pencil (#11 wire) the current would be ©,0.01 amperes. 


Summary 

The actual behavior of a nuclear burst radiation 
source has been modeled by a monoenergetic (1.5 MeV) point 
source of gamma rays. Careful analysis would have to 
consider multiple scattering, actual spectrums, etc. The 
time behavior of the source model is mainly an exponential 
rise and slower exponential decay. If X-rays are of concern 
they are assumed to behave like the gamma source but with 
the energies of 3.8 KeV, the average energy of a l KeV 
blackbody. Neutron time behavior contributes later parts 
of the gamma source model; should a Simple extension of 
neutron behavior be necessary the shell model might be used. 

The gammas and X-rays are treated like a virgin flux 

Eq (63)) from a point source. Of all the photon removal 

processes the Compton effect is the most important means of 
vropelling electrons forward. The electron's average forward 
energy is about half the incident photon energy--graphed in 
Fig. 15. 

Equation (78c) relates the Compton current to the 
photon flux. This equation will be used in the next chapter 


where we consider the symmetric burst in air. 
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IV. Electromagnetic Pulse from 


a Symmetric Burst 


Introduction 

The symmetric burst is most closely approximated by 
subsurface bursts and low altitude bursts that do not inter- 
act with the earth's surface. The field produced is an 
electric field inside the expanding sphere of energetic pho- 
tons from the burst. There will be no magnetic field inside 
ox outside of the expanding sphere and no electric field 
outside that sphere. 

In this chapter an approximate model is used to illus- 
trate the use of Maxwell's equations to find the fields 
inside and outside a sphere of current. The model concludes 
With some discussion of its validity. 

Following this we wili discuss the Compton current and 
conddetton currents. Since the conduction current depends 
on the ionization of the air and the mobility of the charge 
carriers, some discussion of these is included. This electric 
field, which cannot grow indefinitely, is studied to find its 
maximum or Saturated value and the radius from the burst that 
the saturated value can reach. Finally we conclude the 
chapter with an example to indicate the magnitude of the 


effect. 


Spherically Symmetric Model 


By studying a mathematically simple model for a spheri- 


cally symmetric nuclear explosion we can see how Maxweli'ts 
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Fig. 17. Model of Spherically Symmetric Nuclear Burst. 
Characteristics include overall neutrality, positive 
charge evenly distributed over one meter sphere, and 
negative charges evenly distributed in shell with 
expanding outer radius. 
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equations may be used. The model has some special proper- 
ties: A constant number of positive ions remain on G euhere 
of one meter radius and an equai number of electrons are 
propelled radially outward in such a way that the outer 
boundary of the electron distribution moves with velocity y 
and the electrons are evenly distributed between f= 1 
and F=K =Vt. Salient features of the model are pictured 
in Fig. 17. The charge density of the positive ions is 
given a volume charge density by the use of the Dirac delta 
function, the use of which is explained in Refs 44, 45, or 
46. The volume cietes density of the electrons is the total 
charge divided by the volume between [=1 andr = B re 
the moving outer boundary. 

Electric Field by Gauss‘ Law. First let us examine the 


ee ee — 


electric field through the use of Gauss' law in integral 


form: 
I. GE-dS x aay (8) 
&. 
Ss Vv 


The surface bounding the region of interest consists of a 
spherical outer cap at fixed radius [(>])an inner spherical 
cap with radius less than one meter, and a conical surface 
joining the two caps to complete the picture. 

First we will evaluate the left hand side (LHS) of 
Eq (8). Over the inner cap E is zero since any closed 


surface entirely within the one meter radius encloses no 
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charge and fields produced by charges outside r = 1 would 
cancel by Symmetry, 

Next we consider the conical surface. We know that the 
outward normal to the cone (which is the direction of ds) 
is normal to the radius vector forming the cone. By syn- 
metry conditions any component of E must be parallel to 


a. > the radial unit vector. Therefore 


[és 


COME 


(84) 


ni 
W 
‘a7. 
"mM 
{0 > 
Q 
| 
ll 
oO 


which states that no E field penetrates the conical surface. 
All that remains is the outer cap of surface area 
>= ree) where G is an element of solid angle and, by 


differentiation, Cos rdw. The LHS of Eq (8) becomes 
Q) 
(LH: [Eas = [er2ow - Er@w (85) 
fe) 
Now we turn our attention to the right hand side (RHS) 
ef Gauss' integral equation, Eq (8). The RHS can be eval- 


uated by considering three volumes: T€ (0, 4 | ;re(4, r); 


and FOOT (if [> G ). These three integrals are 


Wl 
ae Cs ae ee eee 
(RHS)= Lf [-ae-0r GQwdr + 
06 r 


wr 
i} /-@ 2 als 
Lf Sear drdw eq COJaV (86) 
ome | V 
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where AdV= Csadr = redwdris the cone's volume element. 


The result of the integration is 


(RHs)-L|Qw _ Qu 3 _ ff 
el avr aatteaal © (87) 


The electric field is found from Gauss' integral equation 


formed by equating Eqs (85) and (87): 


SIO. . Oreo) 
E S| anr2 Anrs—1) |r one 


We note that if [ > ff, we have the integrals of 
Eq (86) but contributions from the first two cancel at 
Pos fF, and the third integral is identically zero. The 


conclusion: FE = © outside the charge boundary. 


Let us now consider the magnetic field due to the expanding 


sphere of charge by using another of Maxwell's equations: 
LV GF | = oe D-dS + | J-dS (89) 


If we consider H as the sum of three components in spheri- 


cal coordinates, 
H = Hr art He Ag+ Ho ap (90) 


we can find the components normal to the radius vector by 


considering a spherical cap as shown in Fig. 18. We shall 
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Fig. 18. Illustration of the Cancellation of H by 
Adjacent Incremental Surfaces. 
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discuss the LHS of Eq (89) to find intuitively that the 
tangential components of H are zero and then we shall eval- 
uate the RHS to arrive at the same conclusion. 

We can integrate OH: dl around the circular boundary 
of G in Fig. 18. The value of this integral is hy } where 
Hy is the component of H along £ and #£ is the path 
length around the cap. Hy is constant on this path by 
symmetry. Taking an adjacent cap, a)! » we see in Fig. 18 
that at its point of tangentey Hy =- Ay so the sum is 


zero. Therefore over the entire sphere there is no compon- 


A 


ent of H normal to ar and Eq (90) reduces to H = Pes é 
We shall show later that H,. is zero also. 

“We now proceed to evaluate the RHS of Eq (89) to show 
that the components of H normal to ar are zero. We carry 


out the indicated integration without comment. 


Ww 
DedS = $2 2c (91a) 


, | 
_ wort) 9 =) (ie 


4ar¢+ Nt 
J-dS ‘= (0-7 dve— 2P ay (S2a) 
Ss at 
Vv V 
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1 

ye ew ies 

a [Ss dr Up 2ar cs (92b) 
4 


_ A381)» deal 
4-1 ot ee) 


The two parts of the RHS of Eq (89), Eqs (91b) and (92c), 
are the same but opposite in sign so their sum is zero. 
Since the area of integration was specified only as W , 
the length of the boundary is arbitrary therefore H normal 
to an must again be zero. . 

By choosing another surface it might be possible to 
show that H- = @) also, but we choose to show that 
by using Faraday's law. 

Magnetic Field by Using Faraday's Law. In order to 
show that all the components of H are indeed zero we use 
the electric field, previously found, in Faraday's law, 

VxXE ==, Ly - The electric field, Eq (88), can be 


written 
E. = E(r ap (93) 


By carrying out the curl operation on this field 


aA A = A 
ar fra, rsine 


= _ 1 P 3 d 
VARS rsinetor Je oo (94) 
BE <O- “0 
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we get VxE= 0 (94D) 
or pot - O (95) 


This last equation integrates to show that H is a constant. 
But at t=-O,H =Otherefore H=O for all time whether this 
point in question is inside or outside Lavi . 

Discussion. The model we have just discussed is 
inadequate to model the effects found in the vicinity of 
a nuclear burst. A more realistic model would have to con- 
Sider that the current density increases very rapidly, almost 
a step function, and the radius of the region containing a 
Source current is increasing at the speed of light. Also 
the flux of gamma rays which drives the current is decreasing 
due to geometric (VYr2 }) attenuation and photon removal in 
addition to its late time-dependent decrease. For the time 
dependence of Chapter III the gamma flux would be "on" for 
a time equivalent to ~ 2700 meters and the width about the 
peak is % 300 meters. 

The 300 meter peak pulse region propagates outward 
creating a strong electric field that would be limited in 
Magnitude (called saturated) by the increased conductivity 
of the air, a point to be discussed in the next section. Still 


later the driving mechanisms would not be able to saturate the 


field, and later still no part of the driving pulse would be 
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strong enough to saturate the field so the entire field would 
decay by electron and ion conduction. 

A realistic central region, the isothermal sphere, would 
contain no electric field (just as in our model). This is 
true since the isothermal sphere is a highly ionized, elec- 
trically neutral region. The radius, R, in Fig. 19, of the 
central sphere would grow slowly with time. The current flow 
and electric field would be "on" to some large value just 
outside the isothermal system. 

. The qualitative behavior of the electric field is illus- 
trated in Fig. 19. 

Despite the increased complexity of this process there 

would be no electric field external to the gamma wavefront ' 


and no magnetic field anywhere, just as in our simple model. 


Conduction Currents 

Electromagnetic fields in general and the EMP in 
particular are generated by flowing currents or charge 
distributions. In a nuclear burst the separation of a 
Compton electron and its parent ion creates an electric 


field which will cause currents in a conducting medium. 


The total current density, a is then the sum of this 


conduction current (ge } and the Compton current (Ke } or 


J: J + J. (where we ignore all other form of current- 


polarization, etc.). Using Ohm's law, J = oE > we have 


we 


I oe (96) 


where QO is the conductivity {in mhos/meter). 
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The conductivity depends upon the number of charge 
carriers available and their freedom to move. The initial 
Number of charge carriers is governed by the number of ion 
pairs (an ion and an electron) that the energetic Compton 
electron produces as it is slowed in a medium. The very 
mobile free electrons carry most of the current but they 
Tapidly become attached to atoms which are much less mobile 
than electrons. The conductivity is then the summation of 


all the carriers times their mobilities: 


o(F,t)= YEP LINE tg, (97) 
kel 

where fi is the mobility (velocity per unit field in meters? 
volt~} second); Qi the charge on the jth particle type; 
and N; the particle density of the th species (in particles- 
meter~°). The mobility can be function of the electric field, 
but we shall ignore that dependence in this thesis. 

In order to see the importance of the conductivity to 


any analysis of the EMP let us find a differential equation 


for the electric field in the symmetric burst. In a symmetric 


burst Maxwell's equation VxH = J + a8 becomes 
J +cb+ EAS - 0 (98) 
c at 


since there is no magnetic field. Only the conductivity and 


the driving current J. control this equation. tS terms are 
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interrelated since the Compton current ionizes the air which 
then conducts with a conductivity @ . Fig. 204 illustrates 
the interrelationship of the processes involvea.48 

Because the parameters N, in Eq (97) are net particle 
densities we will study particle production then particle 
attrition. Later we shall consider the mobilities of the 
various particles. 

Charged Particle Production. The charged particles whose 
motion is the conduction current are produced as a high energy 
electron loses its energy through collisions in the medium. 

The process, illustrated in Fig. 20b,begins when a 
photon, shown incident from the left, propells an electron, 
the primary electron, through the medium. *? The primary 
electron with kinetic energy qe loses most of its energy 
in collisions with neutral atoms along its path, some of 
which become ionized. The secondary electrons do not contri- 
bute to the electric field because they do not move far but 
this massively ionized region does form a conduction path 
along which currents flow to neutralize the electric field 
between the primary electron and its parent ion. 

As the primary electron loses energy its velocity 
decreases (which is important in magnetic field interactions). 
The rate of decrease in velocity (or energy) is greatest 
toward the end of its range and at the end of the range the 
rate becomes erratic. This erratic energy loss causes the 
primary electron's range to vary 15 percent for a given 
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initial electron velocity. However, we assume that the 
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Fig. 20. Interrelationship of Field and Currents and 
History of Electrons in Air. In A we show the inter- 
relationship of fields and currents. in B we illustrate 
a photon propelling an electron forward. The subsequent 
electron motion generates an electric field and an 
ionized conduction path. 
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ion pair. production rate and orimary electron velocity are 
constant until the electron abruptly stops at the end of its 
range. 

The mechanism of energy loss is the single ionization 
of an atom along the primary electron's path. This energy 


$1 


loss is quite constant from element to element. We shall 


take this energy loss to be a 32.5 electron-volits/ion-pair 


for an ~ 1 MeV electron in air>*. Therefore the number of 


ion pairs formed by a primary electron is 


Number of _ T_(108) 


Ion Pairs ~— 32.5 (99a) 


3.08 x10" aan On. pales (99b) 


I MeV electron 


where Te iz the electron energy in MeV. The number of ion 
pairs created is also the number of electrons and number of 
positive ions produced. We shail deal with complete ex- 
pressions for the production and loss of charged particles 
after considering the loss mechanism. 

“Air Chemistry. The loss mechanisms are more than the 
neutralization of the charged particies by the recombination 
of an electron with an ion. Under air chemistry we will 
consider the reactions between ions and between electrons and 
ions. Each reaction has a specific rate constant which may 
be affected by temperature, environment (e.g., humidity), 
and mixing that takes place. Although there are many types 
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of reactions, we shall only consider four of them as 
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Table 1 


Hypothetical Parameters for EMP Studies 


Burst Data 


Yield Y = 1 KT* 1 MT* Fissions/Kiloton 
1.5 x 1045 
Prompt X-Rays 0.7 Y 0.5 Y 
X-Ray Temperature 
Prompt Gamma 0.01 Y 0.003 Y 1 Kev 
Prompt Nevtrons 0.01 Y -- Average X-Ray Energy 
3.8 KeY 


Delayed Gamma 0.02 Y 0.03 Y 
Average Gamma Energy 


Delayed Beta 0.02 -- 1.5 MeV 


Debris 0.25 


Growth Constant « 108 sect Decay Constant @ 107 sec”? 


Radiation Data (See also Table 3 for range information) 
Average Compton Electron Energy--0.74 MeV 

Average Forward Energy of Compton Electron--0.62 MeV 

Ion Pairs Formed/0.7 MeV Electron--2.0 x 104 


Environment (See Table 3 also for Air Density) 
Geomagnetic Field--5 x 10-5 webers/m- 

Earth density may be represented by 1060 air density 
Air is 0.7 Nz and 0.3 05; particle density at STP 


2.55 x 1025 m-3 


we mem mw ww we mK He RH ew we wwe ee ee wee 


*Dickinson, H. and P. Tamarkin. "Systems for the Detection 
and Identification of Nuclear Explosions in the Atmosphere 
and in Stace." IEEE Proceedings, 53:1921-1955 (December 
1965). p 1932. 


Other data referenced elsewhere in this thesis. No 
Classified Sources Consulted. 
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important in relatively dense atmospheres; all calculations 
will assume standard temperature and pressure, (STP). 

The four reactions are 

a. The source reaction: A neutral atom receives enough 


energy to become ionized 


X+ AE —+ X+e- (€.g,O> AE = O5 +e-) (100) 


This process was discussed in the preceding section. The 
product of the reaction rate for the production of Compton - 
electrons, Re = ma , » and the number of ion pairs 
produced per primary electron, Eq (99a), is the source tern 


for electrons and positive ions: 


T4108 Te (0%) 


32.5 ah 


S.=S, = PPFD) S 
Inserting numbers from Table 1 into this equation gives 


S.Gs 6:67 x10 20Z00") BS + (102a) 


= 1.5x 10° b (1028) 


where pis the gamma flux. 


b. The electron attachment reaction, 


X+e+M——+>X +M (eg, Oo; +e+M+>+O+M) (103) 
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usually. procedes in the presence of a third body, M > 
because the third body is needed to conserve momentum in 

the reaction between a light electron and the heavy X atom. 
For air the 02 reaction is extremely fast and is the one we 
shall consider exclusively. As an aside we note that the 
reaction rate constant or attachment frequency, xy! > has 


units meter® 


-second™+. Often this is altered by multiplying 
this figure by the number of M molecules present and by the 
number of 05 molecules present (which is equal to a volume 
fraction of M . This is allowable because the fraction of 
0> which becomes ionized is negligible). The attachment 
frequency is uy = 7.4 x 10°39 m®-sec' which at STP in air 

is equivalent to ¥ = % No, Nav or = (0.2N,)N, 

or, finally, Yu. =1.02x10° sed where Nw is the particle 
density of air at STP, 2.68 x 1025 particles-m-3, 54 


c. The electron-ion recombination reaction 
X+e+M—->X+M (eg, NB +e+M— No +M) (104) 


is shown as a three body reaction because of the momentum 
consideration. The reaction has about the same reaction 
rate for Ny and 05 so we will consider the reaction as 
happening for any particle in air. The reaction rate we 
will use is oy = 7.5 x 10738 6 ~ see-} or 

‘of = 2.0 x 10-12 m>-sec7+ where we have absorbed the 


= 
third body term into the reaction. rate constant.>4 


74 


ee eee 


d. The mutual neutralization reaction 


x + ob aK, € x (e.g., O03 +N5> O35 +N.) (105) 


will be assigned the reaction rate ; =23x107 sec" gor all 
reactions whether X and Y are alike or different.°* 

(Other reactions are important in the upper atmosphere 
but unimportant in the denser regions. For a complete 
description of the many reactions see Freyer?> or Whitten 
and Poppoff.>6 Usually the change in importance of a reac- 
tion is a function of the number of reactants, i.e., air 
density. This is especially true where a third body is 
necessary.) | 

The number of each species of charge carrier 1S governed 
by the balance Setween loss and production. By using the 
continuity equation 

(gabe oF change) = (sources) - (losses) - 

(leakage) (106) 
we cen write a differential equation for the number of each 


species.°/ 


We shall assume that the leakage term is zero, 
i.e., no mixing or transport. With this assumption the 
continuity equations for the electrons and positive and 
negative ions are 


aN 
a = Se — KNe —%NLNe oe 
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ora = Se —&.N,Ne — o% N,N_ (108) 
on = YN, — oy N,N_ (109) 


where Se is the source of electrons, Ne is the number of 


electrons/m°; N. is the number of positive ions/m°; and N_ 


be 
is the number of negative ions/m>.>8 As always, the change 
rates are the product of the number of each reactant and 


the reaction rate constant. The fact that charge is conserved 


gives us the equation 
N, = Ne +N_ (110) 


which implies that one of the three linear first order 
differential equations (107), (108), or (109) is not an 
independent equation. 

The three equations were integrated in Appendix B 
after making a simplifying assumption which we look at now. 


Let us alter and insert numbers into Eq (107): 


ave = Se — No(%+ &,N,) ; (11lla} 


= Se - Nol.02x108+2.0x10'-N,) (111d) 
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For <r, N, to be comparable to % , the positive ion 
12 


density must be vl m~>, We state that because the 


positive ion density is below this figure we wili ignore 
XRN, - More justification -for this assumption is pre- 
sented in Appendix B.59 

The integrated electron density equations for various 
time intervals in the source function are 


Nft) = SAD) (ext — ext) ,te(O,T,) (122) 


xT | 
Nt) = Seite | eat Bea yemt-0| steC a) -t222b}* 
Net) = Nett, ) et - 1) a EAT) tea) (126) 
Nt) = NAT) e @Ct-Ts } AST. (i128 
where 


oh 
tT 
wv 
" 


Original source strength (electrons/sec-m°) 


«x, = Time constants of source (sec). 


oO 
on 
rT 


Interval of building source strength. 


vV—cr__ 

io | 
em 
~~ 
1 


= Interval of decaying source strength. 
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2) = Interval of constant source strength. 
t ri = Interval of zero source strength. 
NaC] ) = Electron density at interval's start. 
i, %j,%y- = Reaction rate constants (meters™/sec). 
For times greater than the lower boundary of the 


particular interval the equations become 


N¢t) = Seitle te(T.T) :ToO | aise 


Neate Solr) satel) 


- Ww jte(T, | T) 27, T, (113b) 
5 
Ne = xe ; 
a fr ; ; 
Ary Se te(,T,): GOT G13) 
oy 
Nee 8 17; (1132) 


Rather than discussing the behavior of these equations 
we present the curves in Fig. 21 found by integration of 
Eqs (107) and (108) for early time. The late time behavior 
is estimated from the form of Eqs (113 c, d). There are 
several facts to be noted in these curves: (1) The positive 


ion density has 2 maximum around 1029 ions/m> so the earlier 
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approximation is adequate; (2) The electron and positive 
ion behavior is similar for t< 2 x10 Ssec and (3) The 
electron density is within a factor of 100 of the ion 
density until ~ 107© seconds. This latter point is impor- 
tant as we discuss electron mobility versus ion mobility. 
Mobility. The mobilities of the various charged 
particles are functions of the air density, the applied 
electric field, and particle characteristics, e.g., mass, 
ionization, etc. We shall ignore all variations except the 


difference of mobility between electrons and ions; &9 
Electrons: be ~1.0  meter2-sec-!-voit-} (114) 
Ions: pa, 107 meter2-sec~)-yoir-! (2115) 


By writing the conductivity, Eq (97), so the summation is 


the electron contribution plus all the rest 
T(t) = UNAte *) Li Nig, (116) 
4 
we can see that 


c(t) = WN te (117) 


is a gocd approximation if the electron density is within 


two orders of magnitude of the ion density. 


$0 


reenact 


The balance of this chapter will look at the fields 


around a symmetric burst. The development follows the AFIT 


weapons effects lectures.17 


Electric Field at Early Times 


We are now able to compute the electric fieid around a 
symmetric nuclear burst. We shall look at the very early 
field and a later, nearly constant field. Because the 
initial conductivity and electric field are zero, the early 


Compton current is much greater than the conduction current: 


Je >> cE (118) 


This allows Maxweli's equation Jo+oE+ gD = © to be 


simplified to 


JE 
ee. as (119) 
€ at 
or 
gE = esis cea (120) 
cs YX é, 
Y 
where 
aa mall Cas x 
®, Pee e a, (121) 


The differential equation for the early electric field 


is then 
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= -KY 
dE. é" MERe e pxls 
eae Ne ~ oF vee 


With the boundary conditions EF =O att=O . The soiu- 


tion is . 
E= Eh" MERe _e 1 (xt_4) : Be (123) 
b4Tr* dy SH on of) 
where 


dy Gamma mean free path (meters) 
MIS Electron mean free range (meters) 

«x Gamma e-fold growth rate (sec74) 

© Electron charge (coulomts) 
, Permittivity (1/367 x 10° farads/meter) 

Saturated Electric Field. Shortly after the previous 

equation becomes invalid another assumption can be made which 
leads to a simple equation for the electric field. At some 
point in time the electric field will be unable to grow 
because the conduction current will just oppose and cancel 
the Compton current.©9 at this time the Ee field will not 
vary in time, i.e., 2K =O » and we have what is called the 
saturated electric field. . 


From Eq (98) we have 
ae. - te (124) 
Tg 


The Compton current is given by Eq (68}. The conductivity 


82 


ee 


GNE/PH/71-4 


may be approximated by the electron contribution only since 


100 Ne>N, for t< zp >» (before recombination or attach- 


ment remove the electrons), i-e., 


c FuNe 


(117) 


For Ne we have, by Eq {ll2a), that the number of 


electrons is 


so the conductivity becomes 


o (Ft) = KOOP) 


32.5 49r* (x +7) 


x“ 


eur et _ ext) Ue e 


(112a) 


(125) 


where So the electron source term is given by Eq (103). 


The saturated electric field is then 


= Sr Ree iit 
“OO ane i eo 


Ay §€ 


| Eo? eG ae ‘et thy [Le 
B25 41 r (Xx +2) 


where we have used Eqs (68) and (119) for Jc and O- 


n~ 


a 


4 


| 


(126) 


[a4 
respectively. By clearing fractions and dividing by © 


we get 
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= MFRe 32.5 (Z+a) aA 
ss —- Sa an (1275 
(1 7 e| TeG08) Be 
en* 
Now e&+x)t_.o as time increases since Y ~10°P* 


(where P is pressure in atmospheres) and % x 74 3; there- 


fore the saturated electric field can be written 


MFRe 32.5 (2+%) 4 


(128) 


which is approximately correct for times up to tT: » the 
gamma peak. 61 Notice that Ex, is not a function of the 
weapon yield. 

Maximum Radius of the Saturated Field. In order to 
find the maximum radius of the saturated field we could 
substitute ae for t in Eq (123) since the maximum radius 
occurs when the gamma flux is a maximum. Instead we shall 


substitute the total gamma flux, as » found by integrating 


the flux output fron t = O to t= le ot 


v 


TT 

ica 
Q 
ot 


(1294) 


oo 


a 
[e dt (129d) 
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P, aT 43 
—— < ~1) 


(129c) 


The equations for the electric field, Eqs (125) and 
(128), become the following when the field is saturated to 


its maximum radius: 


= Gy S_, ae aoe E —= Care zal) (130) 


_ MFRe 325(+«) 


(128) 
Te (108 3 | 
e } be 
<7, vA . 
Now we can replace Re =1) ot by “¥" 5 equate the 
(J 
right hand sides; and solve for the terms containing [ ee 
2 6 
p2ek® . YY eTeO") Me Risk 


47 AySo 32.5 (+x) 


This last equation is solvable by iterative procedures. 

Example 3. To illustrate the magnitude of the satu- 
rated electric field and its maximum radius, we will work 
an example for a 1] MT burst at an altitude of 10 kilo- 


63 


meters. Parameters in Table 1 and in Table 2 will be 


used. 
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Summary: 


The symmetric burst has only an electric field E and 
that only within the gamma pulse wavefront. The magnitude 
of the field is interrelated with the air's conductivity and 
the conduction current. The conductivity is the result of 
the primary electron ionizing the air as it is slowed. The 
behavior of the conduction path is affected by the rate of 
removal of charge carriers especially the removal of electrons 
by attachment to oxygen. 

At early times the electron mobility is so much psreater 
than ion mobility that only electrons need be considered. 

By using certain assumptions the early electric field, 
Saturated electric field, and the maximum radius of the 
Saturated field were found. Values for an example were 
w4x10* volts/meter for the autiieared field and its maximum 


radius was v 2 Axi03 meters. 
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V. The Radiated Electromagnetic Pulse 


Because a radiated electromagnetic wave requires both 
E and H fields, the symmetric nuclear burst could not 
radiate a signal (there was no magnetic field produced). 
This is not the case for bursts which are asymmetric or occur 
in the presence of an external magnetic field. These bursts 
generate a radiated EMP through three main mechanisms: the 
asymmetric flow of Compton and conduction currents due to 
differences in the density of the environment; the asymmetric 
flow of currents due to the interaction of moving electrons 
with a magnetic field; and the hydrodynamic exclusion of the 
earth's Magnetic field. As we shall see, a single mechanism 


is usually dominant at a given height of burst. We shail 


study each mechanism separately later in the chapter. 


Height of Burst 

In order to give the reader some insight into which of 
the mechanisms may dominate as the ue ce of the EMP we shall 
qualitatively discuss the variation in EMP generation mechan- 
isms with height of burst. We will divide the range of 
burst heights from surface to far space into five regions. 
In each region we shall discuss the asymmetries present and 
categorize the radiated EMP as strong or negligible. 

Surface Burst. (See Fig. 22-1A) In the surface burst the 
main EMP generating mechanism is the asymmetry in the Compton 


and conduction currents caused by the presence of the earth. 
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These aSymmetric currents radiate much like a dipole antenna. 
There are other asymmetries present, e.g., variation in the 
earth's surface and the presence of the geomagnetic field, 
but their effect is minor compared to the asymmetry due to 
the earth. If we treat the earth as air with 1000 times the 
density of air at STP we find that the mean free range of 
the gammas is 300 meters in air but only 0.3 meters in 
earth. We conclude that there will be a strong EMP from 
the surface burst. 

Air (~3-20 km). (See Fig.22-1B) In the relatively 
dense lower atmosphere (3-20 km) X-rays and gamma rays 
have ranges which are short compared to the scale height of 
the atmosphere (~7 km). This means that the range of each 
wiil not vary significantly with direction. Also because 
of the air's density, the Compton electron's range is too 
short tS permit a Significant transverse current to be 
created by the electron's turning under the influence of the 
geomagnetic field. ®4 

For Example 4 let us find the asymmetric attenuation 
around a nuclear burst in this region at 10 kilometers. We 
Will calculate the e-fold distance for upward and downward 
travel of gamma rays due only to Compton-effect attenuation. 
(The e-fold distance is the distance for a value to change 
by a factor of e or lf/e from an original value, c.f., two- 
fold decrease.} The macroscopic cross section for 1.5 MeV 


3-1 


ammas is 2.5 x 1079 m2-xg7} (= ) or 3.1 x 10° n at 
g & m 
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VERTICAL AIR MASS 


% 


6 


Oo {oO 20 30 40 50 60 rc) 
ALTITUDE (tm) 


Fig. 23. Vertical Air Mass. The mass of air per unit area in 
a vertical column extending upward from an altitude (in Km} 
above sea level. The example is entered at 10 km, proceeds 
through a change in atmospheric depth. of 400 and exits at 

11.1 km (See text). 
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sea level, and the microscopic cross section is 1.2 x 10728 


n*/particle.°> 

By entering Fig. 23 at 10 kilometers and finding the 
altitude change required to cause the vertical air mass to 
change by L/LLm {= 400 kg-m7~7), we can find the e-fold 
distances upward and downward. This procedure is indicated 
for upward attenuation by the arrows on Fig. 23. The e-fold 
distances, ~ 1.1 km upward and «0.9 km downward (found 


by using A(kg-m 25 into denser air}, show in: this example 
that variation in air density causes little asymmetry. 

The minor asymmetries, e€.g., air density variations, 
and geomagnetic field, do cause a radiated signal which can 
be detected over long ranges, but the overall conclusion is 
that the radiated EMP from a burst in this region is weak. ©” 

Air (~20 - 110 km). (See Fig. 22-1C)The air in the 
region from ~20 km to ~110 km is thin enough that the gamma 
rays and neutrons can escape the earth's atmosphere before 
interacting. A space-borne gamma ray detector above a nuclear 
explosion within this region could probably detect ir_68 
The rest of the gamma rays and neutrons interact in the 
atmosphere around 2C€-30 kilometers above the earth. The 
Compton and conduction currents formed at this 20-30 km 
height are initially moving outward from the burst but are 
turned by the geomagnetic field, and these currents, which, 


because of the turning, are now moving transverse to the 


radial direction, generate a radiated EMP. There may be a 
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smali contribution to the radiated EMP by the variation of 
X-ray mean free path with direction but this would be a minor 
source. 

For Example 5 let us calculate the e-fold distance for 
upward, downward, and horizontal travel of the gamma rays 
from a burst at 50 km. By using the same procedure used in 
Example 4 we find the e-fold distance downward is 26 km 
(50 - 24 = 26). When we try to find an e-fold distance 
upward we see that there is not enough atmosphere remaining 
to cause that much attenuation. Horizontally the microscopic 
cross section times the particle density gives a macroscopic 
cross section for 50 km of 2.55 x 10-@ m-1. The reciprocal 
of the macroscopic cross section is the mean free path, 
~ 39 km. These three distances are depicted in the inset 
to Fig. 23, 

In summary, the EMP from this region is quite strong 
due to the transverse currents formed around 20-30 kilometers. 

Near Space ($110 km). (See Fig.22-IIA)The downward 
directed X-rays and gamma rays from a near space nuclear 
explosion ( > 110 km) would impinge uron the earth's atmos- 
phere. The X-rays are absorbed strongly around the 
100-110 km altitude, and the gammas rays are absorbed strongly 
around 20-40 km altitude.°?. Both the photoelectrons and 
Compton electrons which are produced have mean free ranges 


that are long enough to allow the electron to be turned by 


the interaction with the geomagnetic field, For exampie, a 
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0.7 MeV electron at 30 km has a mean free range of 204 
meters compared to a turning radius of only 10 meters. 

Whether the EMP produced by the X~rays absorbed around 
110 km or the EMP from the gammas absorbed around 30 km 
is stronger depends upon the height of the burst: a burst 
near 110 km would probably generate a stronger signal from 
X-rays because of the intensity of the X-ray flux: a burst 
far above 110 km would probably generate a stronger signal 
by the gamma interaction around 30 km. 

In summary, the large transverse currents which resuit. 
from the interaction of the electrons with the magnetic field 
in the near space burst produce a strong, far. ranging EMP. 

Far Space. A nuclear explosion far in space can be , 
expected to produce an EMP due to asymmetries in the case 
which are represented by a shield in Fig. 22-IIB (this EMP 
is called the case signal). The case signal is weak, however, 
largely because the nuclear explosion is energetic enough 
to completely ionize its immediate vicinity and fully 


ae : ; 0 
lonized regions are poor radiators.’ 


The case signal which 
does form results from Compton electrons which are propelled 
outward by the gamma rays being accelerated back to the case 
by the net positive charge which remains at the burst point. 
This acceleration produces a signal. Reference 71 contains 

a complete description of the case signal, including the 


effects of shiclding used in an effort to conceal the 


explosion. 
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Altitude Correction. Because there may be occasion to 


work problems involving other altitudes than those given in 


specific examples, we present altitude correction formulas 


from Karzas and Latter: The mean free path for gamma rays 


and X-rays may be corrected by using 


Ay =Ay oe ¢3xi0° 4 


(132) 
Mey, as ee 
ie to (2 = oy v2 (132d) 
and the mean free range of Compton electrons and photo- 
electrons may be corrected by uSing 
MFR. (MFR,)2-234 (133a) 
: B= 5c FA 
MER ce MPa 2 210 “(EL ) (133d) 


- 


where EF is the photoelectron ener in KeV and ° is 
p rgy Yo 
the ratio of sea level air density to air density at 
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altitude. Table 3 presents basic data about the atmosphere 


for the convenience of the reader. 
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On Investigating the EMP 

Before we discuss the three mechanisms which generate 
the EMP we will look briefly at the methods used for ana- 
lyzing the actual situations and then enumerate a few of the 
major assumptions which we use in this thesis. 

In the literature, the analysis of the various possible 
nuclear burst environments seems to start with an sesunee 
model of the burst environment using idealized concepts 
(e.g., monoenergetic gamma sources from a point, perfect 
sonductons. etc.). Maxwell's equations for this model are 
written in the appropriate coordinate system and often in the 
time retarded form. 

Thes# equations are usually manipulated somewhat and 
then simplified by dropping the small or slowly varying 
terms. Despite all this the equations often must be written 
in terms of integrals which cannot be solved analytically 
or left in differential form. Sotutions can be obtained by 
numerical integration on a computer. 

On occasion analytic solutions are found but the solu- 
tions are often so complex that they disguise the interaction 
between parameters and the consequences of varying them. 
Some modeis lend themselves to an analytic solution but 
beyond certain limits they predict illogical consequences 
(the model doesn't apply). The models we will use are more 
unreal than most of those in the literature but they lead 


to simpler sclutions and do expose the physical principles 
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behind the EMP. If a differential equation is more illumi- 
Mating than its solution might be, we do not hesitate to 
leave it in differential form. 

; Some of the simplifications which we will use are also 
used in the more sophisticated studies on the EMP. One 
assumption is that the burst itself is spherically symmetric, 
i.e., no case Signal. Another is that only one mechanism 
is operating at any one time, e.g., only air density 
asymmetries are present. The final simplification which is 
usually used is the assumption of axial symmetry in the 
environment, i.e., no variation around some axis, hence the 
EMP will Le independent of the azimuthal parameter d ; 
Because of this symmetry there will be only three electro- 
magnetic field components in spherical coordinates: 

E and H or’ EL, Eo » and Hg 


radial’ "1ongitudinal’ azimuthal 


(See Fig. 24). 
THE SURFACE BURST 


The first of the mechanisms for generating an EMP, 
asymmetries in the Compton and conduction currents due to 
variations in the density of the environment, is present in 
all explosions. The asymmetric currents which generate the 
case signal (See Far Space, page 95) are made asymmetric 
by variations in the bomb itself. Because the bombs cannot 
be perfectly symmetric a case signal is always present. The 


asymmetric currents caused by variations in the air's density 
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Fig. 24. Electromagnetic Field Components in an 
Axially Symmetric Environment. 


produce-only a small EMP in the lower atmosphere and are a 
secondary EMP source in the upper atmosphere. Because of 
this minor role we shall ignore variation in air density. 

A major asymmetry in density exists at the earth-air 
interface and is the subject of this section. Here photons 
and Compton electrons have very short ranges in the ground 
compared to their range in the air. This causes a high 
asymmetry in the vertical current flow. 

This air-earth asymmetry occurs even for bursts whose 
fireball does not touch the ground because the gamma's range 
carries them beyond the fireball. However in this thesis 
we shall confine our attention to a model based upon a burst 
on the surface where any photons emitted into the upper 
hemisphere will travel radially in the air and where any 
photons emitted into the lower, geod’: hemisphere will be 
absorbed immediately. In this model we take the conductivity 
of the ground plane to be infinite, i.e., the ground is a 
perfect conductor, and it is considered flat. 

Choosing the surface to be a flat, perfect conductor 
allows us to use simple image electric currents to discuss 
the production of electromagnetic fields. These image 
currents flow, as shown in Fig. 25, antiparallel across the 


surface or continuously through the surface. ’° 


The tangential 
electric field on this surface is therefore zero because the 
antiparallel flows cancel as they must on a perfect 


conductor. If the right hand rule-for determining the 
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Magnetic 
Fieid 
Lines Real Currents 


oO 


i) 


Surface 


Image Currents 


Fig. 25. Image Currents 


direction of a magnetic field is applied to the tangential 
currents above and below the surface, we see that each 

current produces a magnetic field which is directed into 

the paper at the surface. Therefore the surface magnetic 
fields is twice that attributable to the real current. We 
could say that a mathematical tool identified as a magnetic 
current has generated the electric fields just above and relow 
this surface. 

That the electric and magnetic fields play a mathemati- 
cally dual role has been stressed throughout this thesis, 
e.g., writing H= AGB (instead of B= for to emphasize 
the similarity with D = ¢E. We shall use this duality to 


find a magnetic dipole later. 
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Three Phases of EMP Generation for a Surface Burst 

The surface burst can be examined for early, inter- 
mediate and late times. Tne current-production method and 
Simplifying assumptions differ in each phase. We call these 
phases the gamma pulse phase, the saturated field phase and 
the ionic conduction phase 

Gamma Pulse Phase. In the very early moments of a 
nuclear burst an intense, nearly immediate pulse of gamma 
rays is emitted. These create currents that flow just behind 
a wavefront that is moving outward at the speed of light. 
Any electromagnetic signals generated can travel no faster 
than light so no signal is observed until the passage of 
the gamma pulse. To successfully analyze a pulse of this 
Nature would require the full use of Maxwell's equations in 
retarded time. However to picture the events we can assume 
the Compton electrons form a radial electric field ina 
sneli within tie gpmma wavefront and above the earth's sur- 
face (See Fig. 26a). Even the electrons just above the 
surface do not nave enough time to react to the fields 
caused by the image currents just below the surface. In 
this case we will use the electric field given by the early 


time equation from the symmetric field analysis, Eq (123): 


MFRe_e(e ot 


(et 1) 2 
a (123) 
oe Gor me r 


E- doo ee 
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™eo= Angle from Vertical 


-_7 ~~ . 


O05 1.0 
Relative Intensity = 


C of Electric Field 


Fig. 26. Surface Burst. In A, the current and image 
current are moving in opposite directions in a narrow 
‘pulse which is moving outward. In B the_effective 
Magnetic current filaments are depicted. In C we 
show the radiation pattern for this magnetic dipole 
radiator .89 ; 
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Because this pulse is moving at the speed of light we assume 
that we can approximate its effect by an instantaneous 
switching-on of the entire source from the burst point to a 
point at which the source currents are insignificant. We 
restate the previous approximation about ignoring retarded 
time by saying that an external observer would see nothing 
until the production from the entire field reaches hin. 

In order to see that only the ground plane contributes 
to the magnetic field, we find the time rate of change of 
magnetic field, IBy, . We do this by taking the curl of 
the radial electric field, but VXE a3 zero everywhere 
except at the ground plane (This is true because we are taking 
the curl of a radial vector which is a function of [f only, 
except where the current switches in direction}. Since cnly 
the ground plane values of the current contribute to the 
Magnetic field, and hence the radiated signais, let us 
replace the radial electric field by magnetic current loops 
(Fig. 26d) and find the magnitude of the magnetic dipole 
moment which is formed and then, by duality again, return 
to the fields in terms of the Hg vector. 

The magnetic currents are found simply by replacing 


the electric field E by Mg 3 the magnetic surface current 


a 
density:/8 
— 4MFR, e 1 aa eee 
M2 a a 
\s Q 1 €, Ait =a 108, (134) 


nn are 


where we have eiiuestos —l\)by f(t) . Because this 

method depends upon variables which vary sinusoidally in 
time, we shall think of f(t) as one component of the Fourier 
series representation of the electric field's time depen- 
dence. °4 We shall not carry out the Fourier analysis indi- 
cated but this should present no conceptual difficulties 
for, hopefully, the reader is acquainted with Fourier 
analysis and realizes that it is possible. °? In any case 

we re-write Eq (124) to clarify later steps: 


ker 


ae = 
Mg =C ae (135) 


where the C replaces the constants in Eq (134) and the 
time dependence is implied, 

In order to find the total magnetic dipole moment we 
Must integrate the dipole moments of filaments of magnetic 
current. A filament of magnetic current K is CK= Msc 
where is now the radius vector in the ground plane. 83 


The total magnetic dipole moment of the source is then 


KS = |] SdK (136) 


where S is the surface bounded by the magnetic current, 
(Fig. 26b) 
R 


or KS = mt M4 (137 
(w) 
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Ey 
hence KS = ala Cc ra de (138) 
This expression must be integrated from zero to the outer 


boundary of the source region and is 


Reet Ce S|. mcg" 


v7 = 


re = (139) 
but since the value at the outer boundary is negligible we 
set it equal to zevo. This equation is just the electro- 
Magnetic dual of a dipole of electric current. The magnetic 
field of a dipole antenna and its image at the ground plane 
is 


A= ee < as (140) 


which is valid for distances far outside the boundary regions 
(See Fig. 26c).84 The variables in this equation are the 
propagation constant k 5 the angle from the observer to the 
plane's normal @ ; the radian frequency of the particular 
component WJ); the wavelength of the particular component vA; 
the permittivity of free space €, ; and the radial distance 
r, Note that the time dependence is still an implied 
function with a particular frequency. 

| The electric field is found by the relation Peat Hy 
which is applicable for fields far from the source {i.e.; in 


the radiated fiela):&5 
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> -—=~2 Ee (141) 
BK Ar 


where we have inserted the time variation (in exponential 
notation). 

Finally we take the real instantaneous value of this 
equation to get the usual form of the electric field 


equation: 


a tC €&w coswt-kn 
E,= f2 Z, 73 nN ici ae 


sine (142} 


where the\/®2 results from the conversion of root mean square 
amplitudes to instantaneous amplitude. The field we have 
found varies sinusoidally and is propagating outward as 
evidenced by the cos(wt—kr) term. This is the variation 
found in the expression for the radiation far from a dipole 
antenna source. 

In order to estimate the magnitude of the field ve shall 


replace C and insert sample values into this equation: 


7 a %& MER ee w@ cos(wi—kr)sine 
Ee 22.0 a9 “ho ae 2 (143) 


= J? Zz. By MER. 2Q7C cos (wit-kr)sin 8 (144) 
40“ ee Gs 


where 2 = 1207 (intrinsic impedance of free space) 
MFR, = 3 meters 
e = 1.6 x 10719 coulombs 
«x = 108 sec? 
Py = 3.4 x 1025 photons 
c = 3x 10° meters-sec? 


Inserting these numbers we get the following expression: 
ES = 4xigl cos(Wt—kr)sine (145) 
r pa . 


This function is true only for fields far from the source. 
If we choose 10 km and frequencies of 1 Mhz and 10 Mhz we 
get peak values of 4.5 x 10% and 4.5 x 104 volts per meter 
respectively for the two frequencies. This model must break 
down as the frequency grows higher. One discrepancy is 
present because we did not actually compute any Fourier 
coefficients when we conceptually replaced f(t) by its 
Fourier series. 

Saturated Field Phase. In the surface burst, at times 
approximately equivalent to the time at which the electric 
field is saturated, the EMP generation mechanism is different 
from that of the gamma pulse phase. The electric field may 
be nearly constant during this time because the current 
density and conductivity have similar dependence on time 
(See Eqs (126) and (127)). These currents which are 


electronic in character (rather than ionic) flow to a large 
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Fig. 27. Surface Burst with Image Currents. The burst 
is axially symmetric about Z. The lines shown are 
sections through surfaces of rotation. The currents 
intersect the X surface normally. The image currents 
shown dotted) must flow in the direction indicated. 

The inset depicts the current model, Eq (146). 
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extent toward the ground plane as shown in Fig. 27.8® the 
current and image current form a large loop of current 
which creates a transverse magnetic field into the paper in 
Fig. 27. This time varying magnetic field creates a time 
varying electric field also in the transverse direction and 
the mechanism for a radiated signal is formed. 

In line with our discussion of the path of the currents 
in this intermediate time, let us propose a current distri- 


bution 
J = 2 SS" ft) a (146) 


which is shown in the inset to Fig. 27. The radial current 
is greatest along the vertical axis and zero along the 


horizontal axis. The image current given by this function 


does not contradict the behavior specified for image currents. 


We will find the H field that this current generates and 
then discuss the radiated field that it would generate. 


For the saturated-field time regime, Maxwell's law 


VXH = Je TAS becomes Ampere's Law 
VxHoe J (147) 


inside of the current region since the electric field is 
constant in time. By carrying out the indicated curl 


operation we get the following vector differential equation: 


Lili 


UNG Tay tase 


_ its 4 (rsineHs)_ 4 d(r sine H,). : 
rae @ Casey = Sr, Dap (148) 


where we have kept only the Hg component of He oby virtue 

of the behavior of an axially symmetric model. By equating 
A 

the a components of this equation we get a differential 


equation which integrates easily: 


dSiINdH,) . _2d. cose sine f(t) 
© 


(149). 
r 
or 
tm 
sine Hy = 2J, site f(t (150) 
er 
or 


Hs = J, sine f(t) Gets 


r 


This inagnetic field behaves logically for the environ- 
ment we have described. The field is zero at © = 0, i.e., 
along the axis of symmetry and the field's magnitude is 
decreasing with radius. The singularity at [ = 0 is remov- 
able in reality since the net current there is zero. After 
we look at the differential equation for the magnetic field 
external to the source region,we shall see that the specific 


current we chose is unable to provide fields that will match 
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the boundary conditions demanded by a radiated field. In 
order to get those radiated fields we must return co 
Maxwell's equations. 


By using the two curl equations, we shall derive the 


wave equation which the electric fields outside the burst 


must obey. 87 Outside the source region there are ao con- 


duction currents 
(152) 


VxXH =6 


Quigg 
clr 


By taking the curl of this equation and then inserting the 
xE from Faraday's law we get 


Y 


expression for 
(153) 


Vx VX = € oe) 


and 
V(V-H)- VH=-<6 se (154) 


or 
es 
(155) 


Wp abe bl a6 
cd 
1 


E oko” 


where we have used H yey B: VB 20) =1CVH and C= 


This is just the wave equation which has solutions of the 


form 
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H = H(x-ct)+H (xe+ct) (156) 


These solutions are inward and outward propagating waves. 


If the form of the H field is separable into space and time 


functions, the space part of the wave equation becomes 


VH+KH =O  (Helmholz Equation) (157) 


i : 89 
where is the separation constant. 


indicated operation on the H field of Eq (151) we would 


If we carry out the 


find that it could be a solution for Eq (157) only if a) 
was identically zero. That is, the radiatea H field cannot 
meet. the specified boundary condition therefore the boundary 
conditions are invalid. | 

To get a correct solution for this period of the 
surface burst, we would have to use the full set of Maxwell's 
equations and conservation of charge. 

Tonic Conduction Phase. The last phase of the surface 
burst, the ionic conduction phase, begins more than a micro- 


second after the vurst.29 


Low frequency, less intense 
signals are the main features of this phase since ionic 
conduction is slower than electron conduction. The current 
pattern is similar to that of the electronic phase and the 


radiated fieid is somewhat the same. Because of this simi- 


larity we shall not pursue the description further. 


Summary cf the Surface Burst 


The asymmetry of the surface burst means that the first 
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approximation to the radiating mechanism is a dipole radiator. 
The real mechanism could be represented by an expansion of 
the currents as an infinite series of multipole radiators but 
the radiation of higher multiplicity than dipole are less 
important than the dipole contribution. 

Assuming that the earth is a perfect conductor does 
not affect the radiated fields to a great extent since the 
earth is a fair conductor. There would be no penetration 
of the electric field into the earth if the earth were a 
perfect conductor, but because it is not, the fields do 
penetrate the surface. We shall consider electric fields 


in the earth in Chapter VI. 
HIGH ALTITUDE BURST 


The second mechanism for producing asymmetric current 
flow is the interaction of electrons with a magnetic field. 
The asymmetric current from this cause is the major source 


of the EMP from a high altitude burst. 


Electron-Magnetic Field EMP Mechanism 

Any magnetic field will cause an electron to turn in 
an amount proportional to the electron's velocity perpendic- 
ular to the field. This does cause cyclotron radiation 
which is a minor effect in this case. The primary EMP source 
mechanism however is due to the coherent turning of many high 
energy electrons and the coherent flow of conduction current 


in the electric field they form. After a high altitude 


115 


nuclear burst, this coherent motion may be present due to 
the geomagnetic field. 

The geomagnetic field has a value which varies in 
direction and intensity over the earth. For our purpose we 
shall consider a region in the upper atmosphere on the 
magnetic equator with an assigned field strength of 
0.5 x 1074 webers/meter2 (0.5 gauss, egs).>! The Compton 
and photoelectrons interact with this field at all altitudes 
but in the dense lower atmosphere the range of the electron 
is so short that the effect on the EMP is secondary to that 
of other mechanisms. In addition the ranges of the photo- 
electrons are so short that they are absorbed in the region 


92 In the 


ionized by the gamma pulse and cannot radiate. 
upper atmosphere, however, electrons do have long enough 
ranges and lifetimes to interact significantly with the 
geomagnetic field. 

The photon radiation from a high altitude burst inter- 
acts to produce electrons in one or two regions depending 
upon the height of the burst. If the burst is above -v110 
kilometers there are two regions of interaction: 1100-110 km 
where the X-rays interact and \ 20-40 km where the gamma rays 
interact. If the height of burst is between % 20-100 kn, 
only one region of interaction around 20-40 km exists. The 
X-rays are not a Significant EMP source within this region 
for the same reasons this mechanism is not significant in 


the dense lower atmosphere. 
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In order to investigate the reason why a thin layer of 
the atmosphere is the interaction region rather than a dif- 
fice’: broad region we ace consider the absorbtion cross 
section as a function of air density. The photon absorbtion 


functionG(y )is : 


ny) 
. - y 
gly) = Ae Ir (158) 


where we assume a plane wave incident in order to avoid the 
- 4 2: 5 _ . 
geometric ('AAqr™ ) attenuation; -& is the total cross 
section; Ye is the burst height, and A isa proportionality 
9 : 

constant.”> The distance CY, Ye } are measured with ¥y 
positive in the direction of increasing air density (or 
downward) and can be measured from any intermediate level, 
Y, » for this discussion we chose Ys to be the interaction 
altitude. : 

Nee fy) can be represented by different exponential 


functions, each valid over a specific altitude range:°4 


IL Ny) (159a) 


Ly ) 


y 
LL N, e744 (159b) 


(159c) 


" 
= 
D 


ET, aa 


where [lL is the microscopic cross section; Ny ) is the 
particle density aty ; N, is the particle density at Yo» 
+H is the scale height at Y, , and foe is the macroscopic 
cross section at Y, - 


The attenuation function is then 


gly)=Ae % (160a) 


Y, (160d) 


| MA - oe 


(160c) 


aH 
ae 


but in our coordinate system and for the plane wave 


approximation Yo is large and negative therefore the factor 


a7 —O0 


pol on = A/ohe MH 


gy) =A e. 


Now the exponent inQ(y) is small when y<O compared to its 
size when Y > QO. In fact Oy) Overy rapidly because of 
the rapid growth of the term in the denominator of Eq (161) 
or, restating this, most of the photons interact in a shallow 


region.95 
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The behavior of the attenuation function means that the 
density of electrons dislodged by the photon flux will rise 
exponentially in a particular region (but the electron 
density will abruptly drop to zero further down because the 
attenuation removes all the photons). These regions of 
high electron production are thin relative te the distance 
to an observer on earth. 

Gyromagnetic Motion. The electrons dislodged by the 


photon flux experience a force due to the geomagnetic field, 


F = e(V¥ xB) (48a) 


where we have ignored any electric fields (See Fig. 28). 
Since the force is normal to the magnetic field and the 
velocity, the electron motion in a homogeneous magnetic 
field B is circular with a constant velocity (if the small 
cyclotron radiation iosses are disregarded). By equating 
mv, rs 


the centrifugal force, to the field force we can 


find the radius of the circle 


2. 
MMe = ev.B (162a) 
. 
g 
mv 
ree (162b) 
g eB 


where V, is the velocity perpendicular to the magnetic field 
and "g is the radius of the circle. From the relation 


w=Wr we can find the radian frequency of this rotation: 
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Magnetic fieid out of © 
paper 


Fig. 28. Electron Motion in Magnetic Field. In 
a homogeneous constant magnetic field B the 
energetic electron from a photon-electron inter- 
action moves on a circle of radius m=“ "7eB 
The electron will radiate a frequency equal to 
its gyromagnetic frequency, cw/= eB/m? . 
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ie 25 (163) 


These are called the gyromagnetic radius and gyromagnetic 
frequency respectively. 

These equations must be corrected for the relativistic 
increase in the mass of the particle and the resultant | 
equations are 


MmyY . w - £8 (164) ° 


g eB : my 


where Y =/1 —3* and & = ve -* For photoelectrons, 


from 3.8 KeV photons, 4 is 0.12 and for 0.7 MeV Compton 


electrons is. 0.9. The relative sizes of these means 
that relativistic corrections may be ignored for photo- 
electrons but must be considered for Compton electrons. 99 
Electron Motion. Let us return to a more detailed look 
at the motion of the electrons. By supposition the gamma 
rays are produced by a nuclear burst high above the inter- 
action layer. The photons interact in a flat interaction 
layer and produce electrons whose initial velocity is in 
the same direction as the photon momentum. Their subsequent 
motion is a combination of linear motion along a guiding 


axis (horizontal iines at 1, 2, 3 and 4 in Fig. 29) about 


*Relativistic relations: T is particle energy in MeV; 
M,C= 0.511 MeV for electrons; W-—-—,.+ 4 : Bai Ky 


Y= Bt ; yew; mC 
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the same guiding axis? 

The linear motion is due to the component of the initial 
velocity parallel to the magnetic field. The only force 
affecting this motion is the electric attraction of the 
parent ion and any collisions with the atmosphere. The mag- 


Th 


“a sr 
tov,cOosS@ when is/2 . Because of this linear 


% 
nitude of this linear motion varies from zero when G = 


motion, the circular motion appears helical to an external 
observer. 

The circular motion is all we shall consider in our 
example. The circling speed of the electron is V% sin@. 
where Q" is the angle between the magnetic field and the 
gamma momentum. Because V, =V, sin 6" the gyromagnetic 
frequency will be affected through the relativistic 
corrections and the gyromagnetic radius will be altered 
through both the velocity and relativistic terms. 

By considering the motion of the electrons as pictured 
in Fig. 30 we can see the distinctive features of the high 
altitude EMP generating mechanism. The figure depicts a 
gamma pulse five shakes (one shake = 1078 second) wide 
(~15 meters) generating electrons which immediately start 
turning. In the time for these electrons to complete one 
orbit, tens of meters, the gamma pulse has traveled a 
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distance the order of hundreds of meters. The electrons 


in this column have a horizontal velocity component which 


varies with a gyromagnetic frequency W(~w~ 10° radians/second) 
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Fig. 30. Electron Motion in High Altitude Burst. A. A 
series of electrons produced by the leading edge of the 
gamma pulse are shown in their relative positions during 
their first orbit. B. The electric field vector between 
the electron and parent ion is shown in Fig. B. 
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ana whose spacial variation has nodes about three hundred 
meters apart. 

Now because the generating mechanism was 2 photon pulse 
moving with the speed of light, any downward radiation 
generated by the topmost electron will be almost in phase 
with the electrons below it so they contribute to the ampli- 
tude of the signal, i.e., there will be constructive inter- 
ference for a downward moving wave. An upward moving wave 
would at some time in its travels encounter an out of phase 
electron and undergo destructive interference. 

Since the electrons are assumed to work together we 
Shall replace the entire interaction region by a deita func- 
tion current sheet oscillating in the x-z plane. The geom- 
etry of the integration we shall perform is indicated in 
Fig- 31 where we have moved the origin to.an interaction 
height. <A current density. J » Which might be used to find 


a vector potential is 


T(rt) = Ta, Sly) sine” coswt (165) 


where 6 is the angle between the photon momentum and 
Magnetic field and where G) is given by Eq (150) or (261). 
Because the expression for the vector potential would be too 
complicated to integrate if the correct expression for sin@ 


was inserted, we shall use 
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OBSERVER 
(near earth) 


Fig. 31. Geometry for Electron-Magnetic Field Model. 
The photons from the burst interact in the X-Z plane 
forming an oscillating sheet of current density J @,. 
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Je, (2,sin 6 +4.cos e)S(y)sin e'cos wt (166) 


where the factor sin © has been inserted in place of sin 6’. 
This replacement assumes that the burst height is close to 
the interaction plane. 

By replacing the time variation with the real part of 


jot 
€ we can write the current function in phasor notation: 


sr A . ’ ‘ - : ast 

Jt) = Glasin © + 4,cose') §(y) sine e! (167) 
(For a discussion of phascr notation, see Plonsey and Coliin?91 
page 311.) 


The vector potential ACT wD oe for a sinusoidal time 


varying potential is 


= = .  Faran oIkR | 
A(t) = A(T)T(t) = ! Sime dvelat (168) 


ATi} OR 
V 
where k - WE is the propagation constant (meter 7+) 
om vector to source point (meter) 
er vector to field point (meter) 
V' source volume (meter) 


R=f-f vector from source point to field point 
(meter) 
It is re-emphasized that the primes indicate the source 
dimensions, the variables of integration. Since the time 


and space variables are separabie we ‘shall drop ejut until 
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the final expressions are formed. 
Using cylindrical coordinates for the volume integral 
the vector potential becomes 


A(T) = Ag.| [fareees- Jo(a, SIN + a,.cos e') 


o 8 450 


2 lo 2 
a sik WA rey, 
x &(y)sine - (169) 
eG 
which after integrating over iy and replacing 0’ with U’ 


where ae de y* and 2u du! = apts’ becones 
27 oO 
Aly) = ad [= de(a, sine'+ A,cos e) (170) 
—jku' 
xe sine’ 


or 


21 


a ku’ 
Aly) = Fah sin + A.cose)sing delet du! 272) 


“< 


Integrating this we have 
Aiy)= —Jo__14 /@ _ singe’| 4 |sin™e’ 
ON Serie 4-[3 A a) a a 
jeu? 4.17 
«|S | (172) 
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Inserting the limits of integration is straightforward 
/ 
except for U->ow . 
In order to evaluate the limit as U—»oo we nust 
F . k _103 
consider the nature of the propagation constant . 


The propagation constant is, in general, complex: 
-, 4 
ke K jk (173) 


1 x ua 
where k is the intrinsic phase constant and k° is the 
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s 7 
—jKu 
intrinsic attenuation constant. This means ehaee? 


can be written as 


i yt 7# gee 
eJku - eiKU e- Ku (174) 


which is the product of a function bounded by +1 and a 


function that monotonically approaches zero as U’> oo. 


- f 
Therefore in the limit asU oo, e JKU would in fact be 


zero. Other than for this explanation, the complex nature 
of the propagation constant will be ignored. 


The complex vector potential is then 


A(y)= Jeo E 7r.\- ey! (175) 
» 47 = i &p ! k ; 
pee a ee 


The complex magnetic field His found by taking the 


curl of A: 


a 


H= B= & VxA (176a) 


oS Oy a 2 


ae: (176b) 


= Ep ac tas) (176c) 


-_4 jk 
Saioae JKY | (176d) 


Restoring the time variation and taking the real part of the , 


resultant expression gives us the magnetic field: 


H = Re 3. chy et (178) 
= a,—b-coswt- ky) (179) 


The electric field can be found from the expression E,=Z HH, 


where Z is the intrinsic impedance of the medium (1207 for 


105 
free space): 


E = a 301J,cos(wt—ky) (180) 
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We have assumed.in the above that the electrons stay 
in phase with the propagated signal. In fact their velocity 
is less than that of light and therefore the phase differ- 
ences should be considered in a careful analysis. The 
electrons were assumed to have circular motion and no motion 
along the guiding center which is patently false but the 
error should not be serious for our purposes. 

We shall return to this equation in order to discuss 
its validity and compare it to results obtained by Karzas 
and Latter after we discuss the magnitude of the current, 
and the times and frequencies plausible from this mechanism. 

Current Density. In order to estimate the current 
density to use in the model we have proposed, we must discuss 
the equivalence theorem from electromagnetic theory. The 
equivalence theorem asserts that if the fields at a boundary 
are replaced by suitable current sheets, an observer within 
the bounded region would not gob any change in the fields 


around him. }% 


Since in most cases the fields at a boundary 
are unknown, the procedure often used is to guess a plausible 
current. This is the procedure followed below. 

The magnitude of the current density on the boundary 
surface is found from the magnitude of the fields at that 
point. We will use the saturated electric field in the 


interaction layer to find the equivalent current. That 


current is given by 


(181) 
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Fig. 32. Compton Current in Magnetic Field. Only 
electrons (A) starting from within S will pass through 
the square meter shown. Electrons outside of S will not 
pass through this surface; Electron B misses the chosen 
surface; C turns altogether too soon; and D is part of 
another loop systen. 


where Zz. is again the impedance of free space (c.f., Ohms 
law, l= EYp ). 197 We see immediately that our development 
gives the magnitude of the radiated electric field as 
Esa fA 

In order to estimate the saturated electric field we 
can use an approach similar to that used in the symmetric 
burst. However here we must alter the depth of the region 
behind the surface which contributes to current through the 
surface (the MFR, in the planar case}. The turning of the 
electron will reduce the length of the region behind the 
surface as shown in Fig. 32. The effective range is a7 
meters for a gyromagnetic radius of 10 meters. However 
because the electron makes multiple passes through this ere 


we multiply the current production by (electron range/P7fr)- 


g 
The saturated field is then 
Sat ~ ra 
7. G08 ) He 
where n = the number of passes through a region 
Es = the effective production length for current 
io = the gamma production time constant 
Y = the attachment rate (V/ ix at 25 km) 


TI)45 5 = secondary electron production per primary 
electron 
ms re mabe ee he eH His th 
fl. = the mobility (=fe,,., 3! where ie: he 


scale height and Z is the vertical altitude. 


The saturated field found by evaluating Eq (182) is 


_ 45% 7)10°)2.5) 


= (183a) 
“(0.7 x108)1.0 x #78) : 
= 25x103 (volts/meter) (183b) 


Finally we have the magnitude of the radiated signal pre- 
dicted by this model from Eq (194): 2500/4 = 620 volts/meter. 
By using the relativistic correction shown in Karzas and 
Latter's development - (1 -2 y7} - the field strength is 
increased ten fold to 6200 volts/meter. 

Discussion of this result will follow in the next 
section. 

Frequency and Duration. The high altitude gamma EM 
may have the frequencies in it that are in the oscillating 
source currents. Since the expression for the redian 
erequency, W= = PMY » contains the velocity in it there 
will be a spectrum of frequencies in the high altitude EMP. 
A variety of frequencies are shown in Table 4 for a variety 
of electron energies. 

The duration of the electromagnetic pulse from this 
mechanism in a high altitude burst depends on the lifetime 
of the retating synchronous movement of the electrons. 
Using distance = velocity x time we have a means of deter- 
mining the lifetime of the pulse sua it is tabulated in 


Table 5S. 
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Gyromagnetic Frequency vs. Electron 


Electron 
Energy 


3.8 Kev 
0.4 MeV 
0.7 MeV 
1.5 Mev 


2.5 Mev 


*Photon energy from 


Radian 
Frequency 
(rad/sec) 
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2.3 x 10 


1.5 x 


Table 4 


Photon 
Energy* 
Required 


Table 5& 


KeV 
MeV 
MeV 
MeV 


MeV 


Electromagnetic Pulse Duration vs. 


Electron 
Energy 


Kev 
KeV 
10. Kev 
0.7 Mev 
MeV 
2.5 MeV 
2.5 Mev 


** Ranges 


Electron 
Velocity 


(B = v/ce) 


Mean 


Range at { 
Kilometers** 


Free 


(meters) 


25 


360 


2500 


40 


204 


40 


204 
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(100) 
(100) 
(100) 
(20) 
(30) 
(20) 


(30) 


) 


found using Eqs (133a) and (133b) 


Energy 


Gyromag- 
netic 
Radiazs 
(nm) 
4.3 
8.0 
10. 
17. 


25. 


Electron Energy 


Lifetime at 


( } Kilo- 
meters 
(seconds) 
1.4 x 1076 
10. 1076 
42. 107° 
1.5 x 1077 
7.5 x 1077 
1.4 x 1077 
6.9 x 1077 
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Fig. 33. Geometry of Karzas and Latter's Calculations. 
The curved surfaces of the earth, interaction layers, 
and magnetic fields were all assumed flat for their 
calculations. This is acceptable because most of the 
field at the observation point is due to the inter- 
actions between burst and observer. (Adapted from 

Ref 111). 
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Discussion. Karzas and Latter obtain a result which 


we present for the non-relativistic case (photoelectrons) 


and without time retardation: 


FE = E,(sin @ cos ol costwt- kz)- 1a, + 
a 


SING Sin@t- kz) a) sien, (184) 


MI 


6 tp MFReA, (185) 


where the angle © is depicted in Fig. 33. For o=% 
which is the case for the model used in this thesis, their 


equation reduces to 


E = E,sin(wt-kz} a," (186) 


The result arrived at in our model is consistent with their 
resuit since sine and cosine differ only by a phase factor 
and both show the field direction trasverse to the line from 
burst to observer. Both models show the pulse lasting the 
lifetime of the orbiting electrons. !9? The minus one in the 
vernlcos(wt kz) -1] is due to considering the motion along 
the pusdine center parallel to the magnetic field. The 
factor E, was found basically from ef and from procedures 
related to the saturated field calculations of Chapter IV. 


Their expression for Fs is 
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Ex Mmwig (MFRe)(32.5) (187) 
eT tO) 4) 


where a} is 1 for Compton electrons and Yo, WE for photo- 
electrons {(~v0.1)3; Mis electron mass (kg); GW gyromagnetic 
frequency (seem!) 5K 210 /Y electron collision frequency 
(sec71); MFR, is the electron range (m); T1108 32.5 

is the ion pair production figure; and (1 -2 ) is a rela- 


tivistic correction, 119 


For the X-ray interaction around 
100 km Karzas and Latter found an electric field cof 

10 volts/meter and for the gamma interaction around 20-40 
km they found 6 x 104 volts/meter which is at least one 
order of magnitude greater than predicted by our model. 


Apparently we must conclude that our simplified model is of 


limited validity. 
SYMMETRIC BURST IN A MAGNETIC FIELD 


The final mechanism for generating an EMP deals with 
the freezing (see paragraph below) of a magnetic field within 
an expanding ball of highly ionized gasses. To avoid other 
mechanisms we assume that the burst is symmetric except for 


the presence of a uniform magnetic field. 


Hydrodynamic Exclusion of the Earth's Field 
W. J. Karzas and R. Latter!1* describe the electromag- 
netic pulse generated by exclusion of the earth's magnetic 


field from the isothermal blast region, Fig. 34: 
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“Surrounding a nuclear explosion is an intensely 
ionized region in which the air conductivity is 
sufficiently high to freeze in the magnetic field. 

As a result, when the blast wave causes the air to 

expand, the magnetic field is carried along. At 

the same time the magnetic field which is external 

to the ionized region is unable to penetrate the 

ionized region and is pushed outward." 

Hydrodynamic Formulation of the Problem. By following 
Karzas and Latter's method, but filling in some mathematical 
steps, we shail find a differential equation suitable for 
describing this generating mechanism. However the tine 
dependence of the isothermal sphere's geometry and the 
conductivity are too complicated for analytical solution 
of this differential equation so we shall proceed to use a 
simplified model which satisfactorily describes the mechanism. 

The total electric field in a region containing a moving 


conducting medium where the motion involves no separation of 


charge, i.e@., not a plasma, is 


ome -F+VxB (188) 
and Ohm's law, J = re assumes the form 
J =c(E+V¥xB) (189) 
Maxwell's equation, 2S =VXE becomes 
P B V J = D 
SS. Sas Ke 190 
at = ) (190) 


Because the displacement current, oDy, . is small in the 


realm where hydrodynamic approximations are valid Ampere's 
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B field 
vectors 


Fig. 54. Hydrodynamic Exclusion of Earth's Magnetic 
Field. The BO field within the sphere of radius R(t) 
is expanded by the growth of the sphere and the B 
lines outside the sphere are compressed. 
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equation holds: 


= VxB aaa 
TJ = (191) 
Fes 


Therefore the equation governing the behavior of the magnetic 


field in this case is 


(192) 


Isothermal Sphere Model. In order to simplify the 
solution of Eq (192) the expanding isothermal region is 
modeled by an expanding sphere which excludes the external 
magnetic field as it grows in radius. The approximation is 
fairly good as long as the relaxation time, 7 , of the 
magnetic field is quite long, that is 

a zo MRC 1 ) RL) > (293) 
Relaxation time (diffusion time in Jackson, Ref 114) is a 
measure of the time needed for the initial magnetic field 
to decay away; in this case T= TH. RLY where R(t), here 
the isothermai sphere's radius, is a “length characteristic 
of the spacial variation of Bu 115 

The isothermal sphere is nearly at rest, i.e., V -O 
at t= O . Therefore Eq (192) reduces to 

$B. VB. (194) 
st Be 
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and, if O-,coinside the isothermal sphere 


in - Oo (195) 


Wa 


or the magnetic field is frozen. 

The problem has been reduced to one of an expanding 
sphere with the original field frozen in. In the region of 
interest outside the sphere, oO =Oand the displacement 
current is assumed to be zero, 2D _o, therefore 


at 
VxB- O (196) 


and Vi B= O {from Maxwell's equation). Thus we can find the 


solution outside the sphere through the use of a magnetic 


potential © 114 we want the solution of 
Vb,= 9 (197) 
which is 
2 
, > Ay rt. Es F {cos 9) f- O72 31198) 
Azo S 


where Ay and By are constants te be found by matching the 
boundary conditions and the FP COs @) are Legendre polynomials. 
The magnetic field is the negative gradient of the 


potential or 
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B=-—UVad, (199a) 
= P| A 4 | & 
Fred Fees a 
Date = Be) (cose) a. ~ 
£=o r 

ee r 

> ls Se se oeece (198) 

a r de 


The first boundary condition is the value of B at 
infinity: 


dD fe Ay -. a A A 
Blco)= Boa = Bicosea-sine a,) (200) 
The presence of COS Sin this boundary condition 

implies f = 1 since PR €ose)=cose, By retaining only 
f = 1 and equating components from Eqs (200) and Q@990Q we 


get 
A t 1 
a: B,cose = —( A.- Be) P(Cos eS ) 
r CP cf 
A : (201) 
a, —B sine =—(A+ zal = TI tEgos 2 9) 
2 2 P— 00 
Equation (201) becomes 
B.cos © = —A,cos 6 (202a) 
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or 


A.=—B, (202d) 


é 


To find the value for B, we must match the boundary at 
the surface of the isothermal sphere. Inside the isothermal 
sphere the frozen field is carried outward with the expansion 


qilé 


of the sphere so the internal fiel is 


7 
W 


2 
Bees eR) r<R(t) (203) 
2 TRG 


Again equating components as in Eq (201) gives 


| ; | | 
A. R(O) - B(2 ) 
4: BSScose = -(- B= 2) cos 0| (204a) 
R(t) a r=R(t) 
or 
2 3 
B, = a, | Rol) Ray (2040) 


The magnetic field for F>R(t) is therefore 


2 ee oe ig 
B --[-B.—8 [Rot BLE? bose a.—  (205a) 


2 ead 
-F Be 2 0) a a esin @)A, (2050) 
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2 3 
B-B4-—Bh- KO) Riv{cose 4 sine a. (205c) 


In order to get a numerical value for the magnetic field we 
follow the work of Karzas and Latter. 

"Numerical results on the radiation and hydrodynamic 
flow for a sea level nuclear explosion give 


J, 
4 4s 


R(t) = 5e120( oe ies (206) 


for O< tKy% < 0.05 » where t is in seconds and Y in 
kilotons. For OO5 q ty 


Ys 
R(t) =R,.,7 50 Y (297) 


Equation (205) holds until Eq (206) breaks down, after which 
the second and third terms of Eq (295) decay to zero. The 


time constant for this decay is, approximately," 


fo Oo R. (208) 


"where OF is the conductivity of the isothermal sphere 


when R(t ye oe A rough estimate of Os, is” 


rs pag x107 (amp/volt-meter) (209) 


a 
"Consequently, foo % Pete is of the order of a few times 
ts : ; ? : : 
¥ seconds. Since this relaxation time is so long, 
other phenomena, such as fireball rise, and radiative and 


convective cooling, will determine the actual time for the 
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Example 5. Some figures for a 200 Kt explosion are 


given in Table 6. 
Table 6 


Hydrodynamic Exclusion Parameters Due 
to a 200 KT Burst.* 
29.2 + 380 1/3 Radius to 10 PSI 


29.2 meters overpressure - 1600 m** 


290 meters B at 1600 2,6 = 779 


= 0.3 sec ~ 1.003 Boas 


T (decay time) ~ 1200 sec 


*Numbers computed from appropriate equations; **Ref 118. 


From the values in Table 6, we see that the ¢cime to 
maximum sphere radius is relatively long. This means that 
any radiated signal would have a low frequency. That the 
actual time for the magnetic field to relax would not be 7 
is indicated by the large value found here: 7=1200 Sec. 

At a radius equivalent to 10 PSI blast overpressure we see 
that the change in magnetic field is only v 0.003 By Which 
is very small. 

Summary of Hydrodynamic Field Expulsion. The signal 
generated by this mechanism may not be of primary interest 
- atomic defense. However it has been suggested that the 
Signal from this mechanism could be detected even for under- 


ground bursts. While true, it has been pointed out that a 
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Simple gaussian loop which would cancel the earth's magnetic 
field near the burst would be sufficient to eliminate this 


Signal. 


Summary 

The mechanisms for the production of a radiated EMP 
vary. We have studied three methods: the formation of 
asymmetric currents because of differences in the density of 
the surroundings; the formation of asymmetric currents be- 
cause the electrons were given transverse velocity by the 
geomagnetic field; and finally, the third, minor, mechanism 
was the "tweak" of the peonaenecic field by the expanding 
isothermal sphere. Although the basics of each mechanism 
are relatively simple the analysis of actual situations is 
difficult because the mechanisms are present in combination 
with one another. 

Simply stated, the models which approximate the three 
mechanisms are a dipole radiator for the surface burst; a 
radiating current sheet for the high altitude burst; and an 
expanding magnetic field for the hydrodynamic exclusion of 
the geomagnetic field. 

The electric fields strengths found were on the order 
of 104 volts/meter in the first two models. This implies 
that the EMP should be of concern for atomic defense, The 
fact that all modeis formed radiated signals may be important 


for detecting nuclear explosions. 
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VI. Propagation of the Electromagnetic Pulse 


The electromagnetic pulse must be studied as a propa- 
gating signal both for detection purposes and to assess its 
potential for damage. We study propagation of the EMP 
Separate from its generation because the EMP generation is 
often studied using specialized coordinate systems and 
assumptions that do not consider the environment beyond what 
is significant to the explosion. Examples of these environ- 
mental factors include the ionosphere's natural electron 
density; the earth's curvature; and the earth's fiaite 
conductivity. 

We will study three situations: a pulse propagating 
through an ionized layer (such as the ionesphere); a pulse 
propagating throush the neutral region below the ionosphere; 


and, finally, a pulse propagating into the ground. 


Propagation Through an Ionized Region 
The propagation of the EMP through ionized regions is 
studied in order to understand several possible situations. 
from a 
For example, the EMP, low altitude burst might be detectable 
by a satellite but indistinguishable from lightning static 


(spherics) if the EMP's characteristics are altered by its 


passage through the ionized region of the upper atmosphere 


‘(the ionosphere). In another situation, a weak case signal 


from a distant nuclear test might reach the ionosphere but _ 


go undetected on the ground because the ionosphere might 


148 


eee te 


absorb the signal. Another situation requiring the study 

of EMP transmission through an ionized layer occurs when the 
EMP produced at high altitudes (100 km) by X-rays must 
travel through the 1tonized layer around 20-40 km which the 
burst's gamma rays would creates. In this situation the X-ray 
Signal might be completely absorbed. 

The remainder of this section considers EMP propagation 
as it might affect nuclear burst detection by a satellite. 
Because the other examples mentioned would be affected 
Similarly we shall not consider them. 

Two factors that affect the propagation of signals in 
an ionized region are the frequencies in the signal and the 
electron density along the propagation path. These are 
important because an electron density is characterized by a 
minimum frequency below which all signals are attenuated (or 
absorbed if the region is deep enough). We will find this 
minimum freonency by studying the propagation of a2 plane 
wave through an ionized region. The plane wave is described 


by the following equation 


F=E,cos(t —kz) (210) 


or 


E = ReallE,eJ@" -kz)| (211) 


where we have expressed the wave nature? aS a complex quantity 


and where W is the radian frequency of the signal and kK is 
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the propagation constant. The propagation constant is 


K = W/L (212) 


where HM, is the permeability of free space and € is the 
permittivity of the region containing free electrons. !19 


This permittivity is given by 


2 
e= €| T= Ne 


Ne 213 
° we, mt (213) 


where N , © , and M, ere the electron's number-density, 
charge, and mass respectively. If Ne > UEM, the propaga- 
tion constant is imaginary and, in Eq (211), eJKZ becomes 
—kz 


uw 
ec where K is real. This means that the wave is now 


being attenuated exponentially: 


E = Eek ~cos wi (214) 


Typical values for the critical frequency (where 
Ne“ = wée,m, ) range from 3 to 8 Mhz depencing on season, 
time of day, and sunspot activity .120 A typical natural 
electron density profile (Fig. 35) peaks around 100 km and 
again at higher altitudes./*! The critical frequency for 
the 1060 km peak in Fig. 35 is 4 Mhz. 

In Fig. 36 we present a pulse shape recorded 44.6 km 
122 


from a nuclear burst. By looking at the frequency spectrum 


of this pulse, !?5 Fig. 37 we see that most of the signal is 


unlikely to penetrate the ionosphere. However, Karzas and 
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Figure 38. 


Ground and Reflected Electromagnetic Yaves. 


154 


Lek es 


any ruisygs ims 


Latter state that the EMP does have a large enough high 
frequency component that a detectable and identifiable signal 
would penetrate the ionosphere. 124 The spectrum in Fig. 37 
Supports this since the 105 Hz signal is relatively large. 
Karzas and Latter also state that spherics would be largely 
attenuated by the ionized layers because of their low 
frequedcies,*-* and therefore they would be of minor concern 
to the detection problem. 

The pulse from a nuclear burst far in space would be 
affected by comparable attenuation below a comparable critical 
frequency. The electron density following a high altitude 
burst ( %100 km) may be so large that nearly all the electro- 
Magnetic radiation propagating into the ionized region would 


be attenuated (black-out effect). 


Ground and Reflected Wave Propagation 

The second type of propagation we shall consider is the 
propagation of a wave within the region between the ionosphere 
and the earth's surface. An electromagnetic pulse in this 
region will propagate as ground and reflected (or sky) waves. 
A ground wave is an unreflected wave in contrast to the sky 
waves which reflect off the ionospheric layers of the 


atmosphere! ?6 


(See Fig. 38). The signal received at a 
distant station is a combination of these ground and sky 
waves anc must be analyzed by unfolding the signals which 


arrived via various reflection paths if any indication of 


the original signai'ts shape is to be found. Since the 
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Electric Field Amplitude (Volts/meter) 


Time (useconds) 


Fig. 39. Theoretical Propagation of Observed Waveform. The 
observed waveform (Fig. 36) was theoretically propagated as a 
groundwave (left column) and as sky waves. The sun of the 
groundwave and a first reflected wave are shown in the right 
hand column. 
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electronic Signature of a nuclear burst is distinctive, 

the study of pulse propagation may be important to test ban 
treaty policing. After the 1958 high-altitude nuclear 
weapons tests there was a report of radios damaged at a range 


of 2500 miles;!2? 


damage which might have been due to ground 
and reflected waves adding. 

In order to show the effect of propagation on a wave 
shape we present the series of graphs in Fig. 39. These 
graphs show the pulse shape measured at 44.6 km as it would 
look if measured at the distances indicated. The graphs were 
drawn by mathematically propagating the signal around the 
earth as a ground wave and as reflected waves and summing 
the resuits.!78 

The pulse in Fig. 39 is realistic but may not be ex- 
pressed as a simple function so we present a smooth function 
which may prove useful to the reader. Wait (Ref 129) gives 
a representation of the exponential source function et 


which he alters to get a more realistic waveform for a light- 


ning stroke current dipole radiator: 
Ss «le 2 | XK, >Xx (215) 


This equation, for lightning, has a peak around 10-20 x 1076 


seconds. The form for this pulse is shown in Fig. 40. The 
dominant mode (the first wave which "fits" the waveguide 
formed by the carth and the ionosphere) of the propagated 


pulse is shown in Fig. 41 as it would look after traveling 
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Fig. 40. Idealized Source Functions Useful for Computing 
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Fig. 41. Shape of the 2- Asec pulse of Fig. 40 after 
traveling 3009 km. The solid and dotted curves are for 


propagation over earth with infinite and finite conduc-~- 
tivities respectively. 
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3000 km. It is to be noted that the wave is oscillatory 
at this distance and that the curve shown for finite con- 
ductivity has lower and more rounded peaks than the curve for 


a perfect conductor. 


The last signal propagation we shall consider is that of 
the EMP's penetration into the ground. This is an important 
problem since many power and communication links are under- 
ground. Previous studies of the effects of lightning on 
underground systems could provide a specific and detailed 
foundation for studying EMP effects (For example, see Ref 
131). We shall not be as detailed in this paper as the ’ 
cited reference is. If we suppose that the earth's con- 
ductivity is good we can treat the penetration of the electric 
field into earth as a skin effect =robiem. 134 The skin depth, 

S » is a measure of the distance an electric signal will 
penetrate into a conductor before suffering an e-fold 


decrease in intensity. It 1s given by 
&. pre 
§ Wo UL (216) 


where OG is the conductivity in mhos/meter; (Wis the radian 
frequency in radians/second; and {1 is the permeability of 


the medium. Representative values of oO for sea water and 


2 134 


earth are 4.5 and 10° If we 


mhos/meter respectively. 


suppose a 104 yolt/meter pulse with a frequency of 104 Hz, 
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we find that the electric field SO meters into the earth 
(the skin depth) is reduced by e~! but is still 4.4 x 10° 


volts/meter. 


Summary 

The ionized layers can attenuate the EMP significantly 
if the electron density is high enough. After a high alti- 
tude nuclear burst the electron density can be very high and 
Will attenuate most signals incident upon it. The signal 
propagated through the atmosphere is altered significantly 
by waves reflecting off the ionosphere. The EMP will pene- 
trate into the earth a distance wnich can be approximated 


using the skin depth, Eq 216. 
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VII. The Reception of the 


Electromagnetic Pulse 


In general we shall gather into this section examples 
of the ways in which the EMP will interact with a system. 
Again we beg off on analytical or truly accurate solutions 
because of the effort involved in their development and 
solution. We remain focused on simple or approximate models 
intended to illustrate not precisely analyze. 

There are two ways an EMP can interact with a system. 
Both invelve inducing a difference in potential which in turn 
causes a current to flow. This current may represent a false 
Signal in an electronic subsystem or, in extreme cases, it J 
may even destroy the conductor or components it flows through. 


First the electric field will induce voltages, i.e., 


V=-/E-df | (217) 
A 


Secondly the changing magnetic field will induce voltages, 
ging g g 


L.@s 5 
Vi-—— (218) 
We will examine the voltage induced by the EMP in three 
cases: (1) the voltage between two points; (2) the voltage 


induced in a loop; and (3) the voltage on a short antenna. 


Following this we will consider the attenuation of the 
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QO volts on line 


Fig. 42. Electric Field between Two Points. 


In A the cow 
has an electric field between her front and hind legs. 156 
In B, A and B are individually grounded and feel a large 
impressed voltage. 


In C there is only one ground and the 
field in the vicinity has no effect. 


manatee an ane 


impressed signal by shielding with a good conductor and with 


a good conductor having a hole in it. 


The Voltage Between Two Points 

An electric field between two points produces a voltage 
which is commonly accepted as a current source. We shall 
remind the reader of the similarity here to the lightning 
puise and then extend the analysis to a model EMP field. 

In Fig. 42a we see the effect of a lightning stroke 
near a cow. A large potential forms across the ground 
near a lightning stroke; about 30,006 volts/meter for 


“ two meters .!5> 


This potential between the front and rear 
of the cow effectively drives a current through the cow. 

In the case of the EMP, which we assume to be a plane 
wave incident normal to the earth's surface, we find the 
open circuit voltage between A and B of Fig. 42b to he 

£ 
Vs {/E-df = Ef (219) 
} 
If the EMP amplitude is 104 volts/meter and § is 100 meters 
the voltage between the points is 10° volts. 

Now if we look at the grounded system Shown in Fig. 42 
we see that the two grounding points have 108 volts between 
them and if the wire joining the two pieces of equipment has 


7 


a conductivity of 5.8 x 10° mhos/meter and is 2 mm x 2 mm in 


cross section the current flow is 
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Te Wp (220a) 
= yVoA (220b)- 
y 
= 10%58x1074xiO) fog (220¢) 
S 2x 10° (amperes) (220d) 


2, 137 Such a 


(where A is the area of the conductor in m 
current may damage the equipment. 

Of course the solution to this particular problem is 
to build equipment with a ground at one point only, as in 
Fig. 42c.158, this is analogous to having the cow stand on 
one foot during lightning storms. 

This example oversimplifies the solution to the problem 
because the ground cable as shown in Fig. 42¢ has a finite 
conductivity and will therefore develop a potential along 
it which may couple to the active wire. Partial solutions 
are to use a three wire system in which the two signal 


carrying wires are inside a grounded sheath (the third wire). 


Voltages Induced by Time Varying Magnetic Fields 


The second means of coupling the EMP to a circuit is 


by the effect of a time varying magnetic field upon a loop 


164 


GNE/PH/71-4 


within the circuit. We shall in this section derive a 
Simplified formula for the induced voltage and then apply 
this voltage to a time varying pulse of the form given for 
a lightning stroke in the last chapter. Values will be 
asSigned to the amplitude of the EMP field in order to work 
examples which illustrate the magnitude of the induced 
voltages. 

Faraday's law, the third Maxwell equation, is basic 


to the study of voltages due to time varying magnetic fields: 


ean , LBs 


ls a 


We shall concern ourselves only with loops which have fixed 


boundaries so we can rewrite the equation as 


Ve — [28 S (222) 
$ 


and if we further concern ourselves with a plane surface and. 
magnetic fields which are not a function of position on the 


surface S, we can say that 


Vem ses (223) 
Tf B = B, f(t) ,» where 
Hi): te Lent ) > KX (224) 


a) 
oO © 


10 10 
Time (sec) 


Fig. 45. Voltage Induced in a Circuit by a Lightning 
Stroke Model. 
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as in the lightning stroke of Fig. 40, then 


Va-5-B« get _emaut | (225a) . 
OF Oe | é : 
=-§$- Bx xe eas a (225b) 
= $-H, ail ae Si esi (225c) 
Fe ‘ | 


The voltage V is negative for the early parts of the pulse 
If we insert 


since ~ >«& and positive later in the pulse. 


values into this equation we have 


(47 x10- "210 %e-*t —25— “ty (227a) 


10% 
V"S 3567 1207 


ae ity (227b) 


= S(13x104 et -25¢ 
= 5x 10° secm!). For a loop 


where = 2 x 104 sec7! and at 
1 x 1072 nm’, about the size of a 


0.1 x 0.1 meters (S 
hand} the field varies between a negative 3000 volts to a 


positive 90 voits (See Fig. 43). 
Let us return to the example of the two pieces of equip- 
ment separated by 100 meters (Fig. 42¢). Suppose they are 

joined by ordinary household electric-cord (two wires, 
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Fig. 44. 


Short Antenna in an Electric Fiela.14} 
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untwisted, separated by V1 mm of insulation). The area 
between the wires is 0.3 n? and therefore the voltages 
induced will be -30,000 volts to +900 volts. The nose 
effective partial solution to this problem is to use tightly 
twisted wires (VU 18 turns /foot) 139 to reduce the net loop 


size. 


EMP Received on a Short Antenna 
The EMP received on a short antenna, which will be 

described next, is an additional example of an electric field 
inducing a voltage on 2a circuit. We use this example as a 
vehicle to introduce the ideas of equivalent circuits and 
mathematical transformations, specifically the Fourier trans- 
form from the time to frequency domain. The Fourier and 
other transforms are used to relate fields incident upon a 
boundary to fields transmitted through the boundary and so 
are very important for advanced study of EMP transmission. 

Following the work of D. S. Wilson let us consider the 
system depicted in Fig. 44a where we have a dipole antenna 
with a load in the center.}49 By uSing a wire that is short 
compared to any wavelength in the EMP we can make some 
approximations about the effects SE the pulsed field. 
Because this short wire could be some part of a complete 
circuit, e.g., the leads on a carbon resistor where the 
resistor is the load, we see that this is a practical model. 
The equivalent circuit for the short antenna is shown in 


Fig. 44b as a voltage generator V(f )in series with the 
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resistance RA ), capacitance CUt ), and inductance LE ) 
of the dipole and a load Zi % The dipole parameters are ail 
functions of the frequency, f . The load Z, has internal 
capacitance oa resistance Re » and inductance Let - 
The voltage impressed is 
% 
V= |EQ)df = ECCy (228) 


2 
2 
Since the antenna is so short the field is essentially con-. 
stant over the length of the antenna. 

We shall find the natural resonant frequency for the 
antenna system and the frequency spectrum of the driving 
signal in order to illustrate some of the techniques useful 
in EMP analysis. The resonant frequency of a circuit is 
the frequency at which the impedance becomes purely resistive. 


The impedance is 


ae JR? =(w@L— cy (229) 


which is purely resistive when WL= Voc or 
Cf = (230) 


The resonant frequency in the model circuit is 


1 
re ae 1 2 


= Sale ee eee (231) 
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The time dependence of the driving pulse we have 


assigned is 


0 gh 30 
E = (Sot ,t €(O,t,) (232) 
P 
Re gto ky 


In order to find the frequency description of this function 


we must use a Fourier transformation! 4? on it or 


“oo 
V= Eth)f - sfecreer (233) 
<0 
t o 
ig [A eleta . fee JOt at (234a) 
y t, 
0 t, 


_ jet i, tp ; 80 | 
Ef te a 1, fei" al Es BIS ee 
t ~4) jw —~wt 


T jot : tp P 
Xe aoe pr al ~ Ee ell (234c) 
tp jo (j us){ jes) jo 


ee ea (235) 
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A. ANTENNA 


qi 


Fig. 45. Current Induced in Antenna and Loop by Pulse. 
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This voltage as a function of frequency is due to the EMP 
coupling into the circuit in question. If the circuit is 
sensitive to any particular frequency we now have the 
signal's strength and can determine the protection necessary. 
Since this circuit is resonant it will support sustaiied 
oscillations at the resonant frequency (ringing) which damp 


out with the time constant!453 


& Pet Riae, (6ee"!) (25°) 


2(Ly Ly) 

An example of this behavior is shown in Fig. 45a. 
Figure 45b is the current response to a magnetic loop 
antenna. it is to be noted that both circuits oscillate at 
their resonant frequencies and that there is a net current 
flow from the magnetic iocop but not from the electric dipole. 
Naturally this spurious noise could cause difficulties. 

Using Ef to approximate the voltage across the dipole 
is valid only if f< XX for all in the spectrum. This 
would be the case if the circuits under consideration are 
short, the circuits within a radio cabinet, for example. 
If however the behavior of a long antenna were under study 
more careful consideration of the phase of the signals, angle 
of incidence, and other factors of the overall process would 


be in order. 
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Shielding Considerations 

The methods by which an electromagnetic field induces 
currents in a wire, magnetic coupling in a loop and electric 
field potential differences, are straightforward in a very 
Simple circuit. In large pieces of equipment, however, 
there are many complexities present. For example, the various 
circuits scatter any incident fields and the scattered fields 
may be weaker or stronger than the original wave. The 
mutual inductance between wires can couple energy received 
on one wire onto another wire which ray or may not have been 
directly affected by the original signal. The literature 
often uses the generic term Radio Frequency Interference 
(RFI) for such phenomena. 

In order to protect the various circuits from RFI a 
shield is usually put around the equipment. A geod shield 
for radio frequencies would be of thick continuous aluminum 
or copper in order to be a good reflector and to contain any 


currents which were induced. 145 


This however is impractical 
because water lines, input power, antenna couplers, ventila- 
tion holes and access doors all proride entry for unwanted 
Signals. These various unwanted signals interact with 
circuit components where they may be rejected doing no harm, 
or they may cause inappropriate actions within the circuit 
such as changing the state of a flip-flop circuit, or they 
may be strong enough to burn out a circuit. 


The effectiveness of a shield against RFI is measured 


in decibels (db). This unit expresses the ratio between 
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two amounts of power: 46 


S = 10 log a (237) 
2 


Since electrical power is proportional to the electric 


vector amplitude squared we have 


2 
S=10log = =20 log Ee (238) 
2 E, 

where = is the incident field intensity and ES is the 
exiting field intensity. | 

‘We shall look at three factors applicable to attenuation 
through shields: the reflection of signals from the shield; 
the transmission of Signals through shields, and the effects 
of apertures. 

Reflection of Electromagnetic Waves. The reflection of 
an electromagnetic wave from a conducting surface is a 
function of the wave's polarization, angle of incidence, and 
the reflector's properties.: For a linearly polarized wave 
imcident normal to a plane surface the reflection coeffi- 


cient is 
(239) 


where ro and ae are the wave impedances of the two mediums. 


Since there may be material properties and thickness for 


which Z, = Z, , there are possible frequencies for which the 
material is transparent to the electromagnetic signal (e.g., 
optical coatings). However for frequencies of interest in 
the EMP this does not cccur. In order to compute the 
reflectivity of aluminum to waves from free space we use the 
relationship 

2-2 


Alo (240) 


R = = 
2.22 


where Zz. is the impedance of free space (377 ohms} and Zz 

is the impedance of aluminum but Z,) is very smaill--therefore 

R& 1. Because Eq (240) is inadequate to find the attenua- 
i] 


tion by refiection from a shield, we present without deriva- 


tion an equation from C. B. Pearlston: 


6 
R =108.2 +10log SXIO (241) 


Hef 


where G is the conductivity relative to copper 

(Oz. = 5.8 x 107), [lp is the relative permeability 

( eis ~ 1), and f is the frequency in H2.147 For 
example let us find the attenuation a signal would undergo 
by reflection from an aluminum sheet. The conductivity of 
aluminum is 3.54 x 107 mho/meter, Mew 4 » and we choose the 


frequency to be 1 Mhz. Equation (241) is then 


BSH RIO" 5 10° 
= - a SR TO" S 
R=108.2 10109 Faye j 106.1 (242) 


a Oe aR 


The signal penetrating into the aluminum will be 
reduced by 106.1 db. from its initial intensity. 

Transmission of Signals through Shields. If a shield 
is thick enough to have structural rigidity, then higher 
frequency signals will usually be attenuated by that thick- 
ness to a point where the transmitted signal is negligible. }4 
The measure of the rate of attenuation of the field strength 


is given by the skin depth which is the distance for an 


e-fold decrease in intensity. Skin depth is given by 


/ 2 
. a rigor aee (216) 


which for aluminum is 2.7 x 1073, 2.7 x 1074, and g.5 x 1075 
meters for signals of 1 Khz, 100 Khz, and 1 Mhz respectively. 
Note that the depth of penetration is larse for lower 

frequencies (2.7 mm at 1 Khz). Pearlston gives the following 


expression for the attenuation through a shield 
rt ltG 
A= 1.33 x10 t ff Rr (243) 


where t is the shield thickness in meters.149 
The total attenuation by reflection and attenuation is 


given by 
S=A+R (244) 


Apertures. Any aperture in a shield can allow signals 


to leak into the shielded area where they can interact in 


177 


GNE/PH/71-4 


Serrated Spring 
Contact Fingers 


Screened 
Vent Openings 


Remove Finish from 


Contact AreasS——, 


Nassua Sy NCEAGUNtD OO TTLA G Qt Haga cr gety 
i) 
8 
thes wing 


i, 
ett ee ee 


BL 5 


Phone Jack Cover 


Screen Under Dial Mask 


‘ Welded Seam 


| Mave Guide Attenuators Construction 
for Control Shafts 


} 
L__ Gasket 


Conductive Gasket 
Under Cover Plate 


Fig. 46. Features of an Effective Shielded Equipment Case. 
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ways previously described. The study of aperture transmission 
is quite complex even for very simple shapes but there are 
some systematic approaches available. The first approach is 
to model the aperture by a radiator of the same shape and 
size, for example, 2 half wave length slot may be modeled by 
a half wave antenna. Another approach is to replace the 
electric field across the aperture by the equivalent magnetic 
current sheet and then find that current's radiation pattern. 
Another possible approach when the aperture opens into a 
cavity is to treat the system as the aperture driving a 
resonant cavity. 

We shall make a few generalized statements in lieu of 
an analytical attack:150 

Any holes into the protected area must be kept small 
in comparison to the operating wavelength. 

For 2a given distance from a single hole the leakage 
intensity is proportional to the cube of the hole dimensions, 
e.g., doubling the diameter of a round hole will cube the 
intensity of the leakage field. 

Wave guide attenuators are useful for large holes, since 
they pass frequencies above a cutoff frequency and attenuate 
all frequencies below that frequency. 

In Fig. 46 we see the features of an effectively shielded 


equipment case. 
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Summary 
The two basic mechanisms for coupling the EMP to a 
physical system are most clearly stated by the equations 


! 
Vz Edf ; Vege 


The EMP is attenuated by shielding which must be thick and 
continuous to be completely effective. Since shielding 
must always have entry and egress holes, engineers must 


design shields and circuits keeping the EMP in mind. 


Final Comment 


This thesis provides a basic description of the EMP 


from source to coupling with circuits. Improvements in 


future work might include finding a Simple, self-consistent 


surface-burst current distribution that will match the 


boundaries conditions imposed by the radiated field equations. 


Another improvement would be the addition of a Simple aper- 


ture leakage calculation to illustrate the physics of the 
problem. For this last problem it is suggest that Ref 152 


be found. 
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Gauss‘ Law 


No Free Magnetic Poles 
Faraday's Law 


Ampere's Law (with 
displacement current) 


Conversion Factors 


12 calories 


KT = 10 

= 2.6 x 1025 Mev 

= 4,184 x 1012 Joule 
Cal = 4.184 Joule 

= 3.086 ft-lb 

= 2,61 x 1073 Mev 

= 3,966 x 1075 BTU! 
eV = 1.6 x 10719 Joule 

= 11610°K 


kb (Kilobar) = 108 newtons/m* = 986.92 atmospheres 
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Appendix B 


Integration of Electron Density Equations 


The differential equation for the electron density, 
Eq (107), can be integrated for any time dependence if the 
recombination term is neglected and if the electron source 
term is separable, i.e., 
SC. t)> Sri) 
Karzas and Latter do this in their treatment of air 


1 Baran also follows this method,* We shali 


conductivity. 
integrate the differential equation for the electron density 
using four forms of the time dependence as pictured in 
Fig. 6. Fellowing the integration we present a brief justi- 
fication for neglecting the recombination tern. 

In this Appendix the notation 0 ( ~~ ) means the opera- 
tion in the brackets was used to proceed from one step to 
the next, e.g., Eq A; 0 (x 2); Eq B means that Eq B was 
obtained by multiplying Eq A by 2. 

From Eq (107) we have the basic differential equation 


(neglecting recombination) 


dn. _ + 

GN. . YN, = S(r)Tt) (81) 
ocxe#t dt ) This makes the LHS a perfect differential. 

d(Nee™*) = s(r)T(tyet* at (B2) 
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O({Integrate--we will assign lower limits later.) 


ner! t 
fa Ne“) = S(r) [roette (BB) 


init. initial 


For the time interval fromt=O tot= 7, the time 


ax ce 


dependence is I= therefore Eq (B3) becomes O(Let 


T(t)s et and apply limits) 
N. € 
Nee*! = suo fe e% at (B4) 


O (integrate; insert limits) 


Noe%t - Str) (eer _4, (BS) 


K+ YB 


which becomes 0 (x et ) 


Ne ete fet est -et (B6) 


which is Eq (112c). 


oC t>Od ;e-’t_,o ) 


Ne = ehh) eo (B7) 


This is Eq (il5a). 
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t€(T, 7) 
In the next time interval, t 2s Se » the time 
Td 
1 


dependence is T(t) = eAlt a therefore Eq (B3) becomes 


Q(Replace T(t) ; write limits) 
ae t 
Id(Nee**) = Sirjeth [na —Th) 6%t at (B8) 

a tL 


= 


4 


where T= S(ryel™%*O 1) 


and where S(ryex | replaces S0r) in order to match 
boundaries and the lower Limit of the LHS is the value of 
Ne el, from Eq (AS). O({Integrate) 
N et S&r) [ea 
€ K+ UY 
«lt Al, 
_ S(rje te ! Ce ~4yt ae | oe 
Z- 4 
QO(transpose; multiply all terms out} 
BA - 4 
ut Sgrje" eh st Qxt 
Nee pc 
a4 
a A, wed, (+ % )T 
Sole e€ e " S(r)e 1 Sk r) (B10 
yA wrx Brn 
_yt T, 
o(xe % ; factor S.(rje~? out) 
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+ 


No sings |S Se 
er u-B Xa 


' e ‘Q. 


a x+y 


elt -7,) eo xy —ZYT | : 
SS (Bil) 


O(Drop last term as very small and decreasing; factor out 


Vv 


/4-G@ ; Clear fractions in last term.) 


«I, 
y,- Sie io “0 _(s08 out -7, | (B12) 
—~3 


% 17% | 


which is Eq (112b). 
O(Let t>>T, ; since Y>> Athe last term in bracket drops 


out leaying Eq 113b)) 


te Cb Ts) 
For the third integration over the interval te(T,, T, ) 


the source is constant or {T(t)= 1 and Eq (B3) becomes 


Nee" t ‘ 
+ sae 
Ne = Sirye%h ooh 0 fox dt (B23} 
UT. 
NiQe* 2 Tt 
where Si m) exh ek ied replaces Sat) 


CCIntegrate; apply limits) 


ean ae —T; ) 
- 1 3 t , 
N,e*" ~NT,)e*® « Sir) ee Zee 


a 


(B14) 
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0(x et 3; transpose; factor out term) 
xT, ~aAT, —T, ) 
Ne Sire ares 
% 
(4 rT.) Lyxt=T.)> 7 
x j-e% ~ 4a ~ N (T2)e 


B15 
S(r)esh Atk u) | 2) 
The factor NT.) could be replaced by the final value from 
Eq (B12) but this adds nothing to the analysis. 

O(Let t>>T, to get the steady state condition with a 


constant source.) 


Niet 
% % 


e (B16) 


which is Eq (113c). 


ed eae 
For times greater than Tz there is no source,! =O 
and we return to Eq (107) in order to reformulate the 


differential equation: 


|? 


> —2N, — &<,N,N_ (B17) 


OQ 
ct 


Again we shall ignore recombination o(xdtyy 9 


=-— pdt (B18) 
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cae rieme oT 


O({Integrate with limits from t = T, tet 


3 (B19) 


which is Eq (112d). For t>>T,we note that N,7Oas we 


would expect with no source. 


On the Validity of Dropping 
Foliowing Barron we find additional support for the 


validity of dropping N,in comparison to 2 by integrating 


the positive ion production from t = 0 to t = Ty assuming na 
losses: 
a 
xt 
N. = S.C re dt (B20) 


(B21) 
Inserting the source strength used to compute Fig. 21 gives 


25 
N = SAX10 — = {e3 - 1] (B22) 
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7 x19 !8 Pest eaye sens (B23) 
meter 


He 


which is below the approximate-validity-limit of 1019, 
For a discussion of the effect of neglecting recombi- 


nation see reference one page 
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